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Abstract 

We consider a class of backward stochastic differential equations (BSDEs) driven 
by Brownian motion and Poisson random measure, and subject to constraints on the 
jump component. We prove the existence and uniqueness of the minimal solution for 
the BSDEs by using a penalization approach. Moreover, we show that under mild 
conditions the minimal solutions to these constrained BSDEs can be characterized as 
the unique viscosity solution of quasi- variational inequalities (QVIs), which leads to a 
probabilistic representation for solutions to QVIs. Such a representation in particular 
gives a new stochastic formula for value functions of a class of impulse control problems. 
As a direct consequence we obtain a numerical scheme for the solution of such QVIs 
via the simulation of the penalized BSDEs. 
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1 Introduction and summary 



Consider a parabolic quasi-variational inequality (QVI for short) of the following form: 



mm 



^-Cv-f,v-nv\ = 0, on [0,r)xM'', ^(r,-) = g on M^, (1.1) 



dt 

where C is the second order local operator 



Cv{t,x) = {b{x),D^v{t,x)) + ^tT{aa''{x)Dlv{t,x)) (1.2) 
and TC is the nonlocal operator 

7iv{t,x) = sup [z;(t, x + 7(x, e)) + c(x, e)]. (1-3) 

In the above, DxV and D^v are the partial gradient and the Hessian matrix of v with respect 
to its second variable x, respectively; stands for the transpose; (•, •) denotes the scalar 
product in W^; S"^ is the set of all symmetric d x d matrices; and E is some compact subset 
of Ri. 

It is well-known (see, e.g., [3]) that the QVI (jl.ip is the dynamic programming equation 
associated to the impulse control problems whose value function is defined by: 

v{t,x) = sup B\g{X'/n+ r f{Xl^^nds+ c(X*L^'", • (1-4) 



t<Ti<T 



More precisely, given a filtered probability space (0, J'-', P, F) where F = {J-t}t^ we define 
an impulse control a as a double sequence {Ti,(,i)i in which {xj} is an increasing sequence 
of F-stopping times, and each is an ^^-^ -measurable random variable taking values in E. 
For each impulse control a = {Ti,^i)i, the controlled dynamics starting from x at time t, 
denoted by X*'^'", is a cadlag process satisfying the following SDE: 

where is a d-dimensional F-Brownian motion. In other words, the controlled process 
^t,x,a gYQiygg according to a diffusion process between two successive intervention times ti 
and Tj+i, and at each decided intervention time Tj, the process jumps with size AX*^^'" := 

i i 

We note that the impulse control problem (jl.4|) may be viewed as a sequence of optimal 
stopping problems combined with jumps in state due to impulse values. Moreover, the QVI 
(jl.ip is the infinitesimal derivation of the dynamic programming principle, which means 
that at each time, the controller may decide either to do nothing and let the state process 
diffuse, or to make an intervention on the system via some impulse value. The former is 
characterized by the linear PDE in (jl.ip . while the latter is expressed by the obstacle (or 
reflected) part in (jl.ip . From the theoretical and numerical point of view, the main difficulty 
of the QVI (jl.ip lies in that the obstacle contains the solution itself, and it is nonlocal (see 
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(jl.3p ) due to the jumps induced by the impulse control. These features make the classical 
approach of numerically solving such impulse control problems particular challenging. 

An alternative method to attack the QVI (jl.ip is to find the probabilistic representation 
of the solution using the Backward Stochastic Differential Equations (BSDEs), namely the 
so-called nonlinear Feynman-Kac formula. One can then hope to use such a representation 
to derive a direct numerical procedure for the solution of QVIs, whence the impulse control 
problems. The idea is the following. We consider a Poisson random measure fi{dt, de) on 

X E associated to a marked point process {TiXi)i- Assume that /x is independent of 
W and has intensity X{de)dt, where A is a finite measure on E. Consider a (uncontrolled) 
jump-diffusion process 

X, = Xo+ r h{Xu)du + r a{Xu)dWu + V 7(^T.- ' C*)- (1-6) 
Jo Jo ^^<, 

Assume that u is a "smooth" solution to (jl.ip . and define Yt = v{t,Xt). Then, by Ito's 
formula we have 



Yt = giXr) + r fiX,)ds + Kt - Kt - C {Z.,, 
Jt Jt 



dW, 



^ I iUsie)-c{X,^,e))fiids,de), (1.7) 
Je 

where Zt = a' {Xt-)DMt, Xt-), Ut{e) = v{t, Xt- + -/{Xt- ,e)) - v{t, Xt-) + c{Xt- , e), and 
Kt = Jq{— — — Cv — f){s, Xs)ds. Since v satisfies (jl.ip . we see that K is a continuous 
(hence predictable), nondecreasing process, and U satisfies the constraint: 

-Ut{e) > 0, (1.8) 

The idea is then to view (|1.7p and (jl.Sp as a BSDE with jump constraints, and we expect to 
retrieve v{t,Xt) by solving the "minimal" solution (Y, Z,U, K) to this constrained BSDE. 

We can also look at the BSDE above slightly differently. Let us denote dKf = dKt — 
J^Us{e)fi{dt,de), t > 0. Then K is still a nondecreasing process, and the equation (II. 7p 
can now be rewritten as 

Yt = g{XT)+ r f{Xs)ds+ r [ c{X,^,e)fi{ds,de)- f {Zs,dWs) + Kt - Kf (1.9) 
Jt Jt Je Jt 

We shall prove that v{t, Xt) can also be retrieved by looking at the minimal solution 
(y, Z,K) to this BSDE. In fact, the following relation holds (assuming t = 0): 

v{0,Xo) = mi{y£R :3Z, (1.10) 

y+ r {Z,,dWs)>g{XT)+ r f{X,)ds+ r / c(X,-,e)Mds,de)}. 
Jo Jo Jo Je 

We should mention that (jl.lOp also has a financial interpretation. That is, v{0,x) is 
the minimal capital allowing to superhedge the payoff IIt{X) = g{XT) + /q f{Xs)ds + 

c{X g- , e) iJ,{ds , de) by trading only the asset W. Here, the market is obviously incom- 
plete, since the jump part of the underlying asset X is not hedgeable. 
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Inspired by the above discussion, we now introduce the following general BSDE: 



Yt = g{XT) + f{Xs,Ys,Zs)ds + KT-Kt (1.11) 

- r{Zs,dW,)- r [ {Us{e)-c{X,-,Y,-,Zs,e))fi{ds,de), 0<t<T, 
Jt Jt Je 

with constraints on the jump component in the form: 

h{Ut{e)) > 0, Ve E -E, < t < T, (1.12) 

where /i is a given nonincreasing function. The solution to the BSDE is a quadruple 
(y, Z, [/, K) where, besides the usual the component (y, Z, [/), the fourth component K is 
a nondecr easing, cadlag, adapted process, null at zero, which makes the constraint (jl.l2p 
possible. We note that without the constraint (jl.l2p . the BSDE with K = Q was studied 
by Tang and Li [21] and Barles, Buckdahn and Pardoux [2]. However, with the presence 
of the constraint, we may not have the uniqueness of the solution. We thus look only 
for the minimal solution (y, Z, [/, K)^ in the sense that for any other solution (y, Z, [/, K) 
satisfying (jl.lip - (jl.l2p . it must hold that Y <Y . Clearly, this BSDE is a generalized version 
of (|1.7p - (|1.8p . where the functions / and c are independent of y and z, and h{u) = —u. 

We can also consider the counterpart of ()1.9p . namely finding the minimal solution 
(y, Z, K) of the BSDE: 

Yt = g{XT) + r f{X„Y„Z,)ds + r [ c{X,-,Y,-,Zs,e)fi{ds,de) 
Jt Jt Je 

{Z,,dWs)+KT-Kt, 0<t<T. (1.13) 

It is then conceivable, as we shall prove, that this problem is a special case of (|l.lip - ()1.12p 
with h{u) = —u. 

It is worth noting that if the generator / and the cost function c do not depend on y, z, 
which we refer to as the impulse control case, the existence of a minimal solution to the 
constrained BSDEs ()1.7p - (jl.8p may be directly obtained by supermartingale decomposition 
method in the spirit of El Karoui and Quenez [11] for the dual representation of the super- 
replication cost of IIt{X). In fact, the results could be extended easily to the case where 
/ is linear in z, via a simple application of the Girsanov transformation. In our general 
case, however, we shall follow a penalization method, as was done in El Karoui et al. |10j . 
Namely, we construct a suitable sequence (y", Z", [/", i('") of BSDEs with jumps, and 
prove that it converges to the minimal solution that we are looking for, by using a weak 
compactness argument of Peng [18] . 

Our next task of this paper is to relate the minimal solution to the BSDE with con- 
strained jumps to the viscosity solutions to the following general QVI: 

dv 

min^- — - Cv - f{-,v, a' D^v) , h{nv -v)^ = 0, (1.14) 
where TC is the nonlocal semilinear operator 

Ti.v{t,x) = sup[v(t,x + ^{x,e)) + c{x,v{t,x),a'' {x)Dxv{t,x),e)]. 

eeE 
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Under suitable assumptions, we shall also prove the uniqueness of the viscosity solution, 
leading to a new probabilistic representation for this parabolic QVI. 

We should point out that the BSDEs constraints have been studied by many authors. 
For example. El Karoui et al. [10] studied the reflected BSDEs, in which the component Y is 
forced to stay above a given obstacle; Cvitanic, Karatzas and Soner [S], and later Buckdahn 
and Hu [B] considered the case where the constraints are imposed on the component Z. 
Recently Peng [18] (see also |19] ) studied the the general case where constraints are given on 
both Y and Z, which relates these constrained BSDEs to variational inequalities. The main 
feature of this work is to consider constraints on the jump component (U) of the solution, 
and to relate these jump-constrained BSDEs to quasi-variational inequalities. On the other 
hand, the classical approach in the theory and numerical approximation of impulse control 
problems and QVIs is to consider them as obstacle problems and iterated optimal stopping 
problems. However, our penalization procedure for jump-constrained BSDEs suggests a 
non-iterative approximation scheme for QVIs, which, to our best knowledge, is new. 

The rest of the paper is organized as follows: In Section 2 we give a detailed formulation 
of BSDEs with constrained jumps, and show how it includes problem (|1.13|) as special 
case. Moreover, in the special case of impulse control, we directly construct and show the 
existence of a minimal solution. In Section 3 we develop the penalization approach for 
studying the existence of a minimal solution to our constrained BSDE for general /, c, and 
h. We show in Section 4 that the minimal solution to this constrained BSDE provides a 
probabilistic representation for the unique viscosity solution to a parabolic QVI. In Section 
5 we discuss numerical issues for approximating QVIs by a penalization procedure. Finally, 
in Section 6 we provide some examples of sufficient conditions under which our general 
assumptions are satisfied. 

2 BSDEs with constrained jumps 
2.1 General formulation 

Throughout this paper we assume that (0, P) is a complete probability space on which 
are defined a d-dimensional standard Brownian motion W = (Wt)t>Q, and a Poisson random 
measure ;U on x E, where E is a compact set of M'', endowed with its Borel field £. 
We assume that the Poisson random measure /u is independent of W, and has the intensity 
measure X{de)dt for some finite measure A on {E,£). We shall often assume that the 
support of A is the whole space E, i.e. 

Ve G E, 30 open neighborhood of e, X{0) > 0. 

We set fl{dt, de) = ^{dt, de) — X{de)dt, the compensated measure associated to n; and denote 
by F = {J-'t)t>o the augmentation of the natural filtration generated by W and fi, and by 
V the cT-algebra of predictable subsets of x [0, T]. 

Given Lipschitz functions b : ^ W^, a : R'^ ^ W^^'^, and a measurable map 7 : 
M.'^ X E ^ M.'^, satisfying for some positive constants C and k^, 

sup 17(2;, e)| < C, and sup |7(x, e) — 7(x', e)| < k^\x — x'\, x,x' ^W^, 

eG-E eeE 
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we consider the forward SDE: 



dXs = b{Xs)ds + a{Xs)dWs + [ -f{X,- ,e)n{ds,de). (2.1) 

Je 

Existence and uniqueness of ()2.ip given an initial condition Xq G M.'^, is well-known under 
the above assumptions, and for any < T < cxd, we have the standard estimate 



E 



sup \Xt\ 

0<t<T 



< oo. 



(2.2) 



In what follows we fix a finite time duration [0,T]. Let us introduce some additional 
notations. We denote by 

• 5^ the set of real- valued cadlag adapted processes Y = (yi)o<t<T such that 2 •= 



E 



SUPo<t<Tl^t| 



< 00. 



LP(0, T), p > 1, the set of real- valued processes {(j)t)o<t<T such that E \ 4't\'''dt 



< 



00; and Lp(0,T) is the subset of LP(0,T) consisting of adapted processes. 
• LP(W), p > 1, the set of M'^-valued P-measurable processes Z = {Zt)o<t<T such that 



\Z\ 



LP(W) 



E 



Io\Zt\'dt 



< 00. 



LP(/2), p > I, the set of V .^-measurable maps U : 0, x [0,T] x E 



such that 



\U\ 



nlo lE\Ut{e)\PXide)dt 



< 00. 



• the closed subset of consisting of nondecreasing processes K = {Kt)o<t<T with 
Ko = 0. 

We are given four objects: (i) a terminal function, which is a measurable function 



9 ■ 



satisfying a growth linear condition 

sup II 



< 00, 



(ii) a generator function /, which is a measurable function / : M" x M x 
a growth linear condition 

\f{x,y,z)\ 



sup 



{x,y,z)GR'^xRxR'i 1 + f| -|- \y\ + \z\ 



< OO, 



(2.3) 
satisfying 

(2.4) 



and an uniform Lipschitz condition on {y,z), i.e. there exists a constant kf such that for 

ah X G R'^, y, y' G R, z, z' G W^, 



\f{x,y,z)- f{x,y',z')\ < kf{\y-y'\ + \z-z'\), 

X E 



(iii) a cost function, which is a measurable function c : I 
growth linear condition 

\c{x,y,z,e)\ 



X M X 



sup 



{x,y,z,e)&.'ixRx9.<ixE 1 + F| + \y\ + \z\ 



< OO, 



(2.5) 
satisfying a 

(2.6) 
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and an uniform Lipschitz condition on {y,z), i.e. there exists a constant kc such that for 
all X e R'^, y, y' G M, z, z' E M"^, e G E, 



\c{x,y,z,e) - c{x,y',z',e)\ < kci\y - y'\ + \z - z'\), (2.7) 

(iv) a constraint function, which is a measurable map h : R x E ^ s.t for all e £ E, 

u I — > h{u, e) is nonincreasing, (2.8) 

and satisfying a Lipschitz condition on u i.e. there exists a constant kh such that for all 
u,u' GR, e £ E, 

\h{u,e) - h{u',e)\ < kh\u-u'\. (2.9) 

Let us now introduce our BSDE with constrained jumps: find a quadruple {Y, Z, U, K) 
e X L2(W) X L2(/2) X satisfying 

Yt = g{XT) + j'^ f{Xs.Ys,Zs)ds + KT-Kt (2.10) 
{Zs,dWs)- [ [ iUs{e)-c{X,-,Y,-,Zs,e))fi{ds,de), 0<t<T,a.s. 

Jt JE 
with 

h{Ut{e),e) > 0, dP ^ dt ® X{de) a.e. (2.11) 

and such that for any other quadruple {Y , Z, U , K) E 5^ x L^(W) x L^(/i) x A'^ satisfying 
(f230]) - (f2TT]) . we have 

Yt < Yt, 0<t<T, a.s. 

We say that Y is the minimal solution to (j2.1(jp - (j2.11|) . In the formulation of Peng, one 
may sometimes say that Y is the smallest supersolution to (|2.1U|) - (|2.11|) . We shall also say 
that (y, Z, U, K) is a minimal solution to (|2.10p - (j2.1ip . and we discuss later the uniqueness 
of such quadruple. 

Remark 2.1 Since we are originally motivated by probabilistic representation of QVI's, 
we put the BSDE with constrained jumps in a Markovian framework. But all the results 
of Section 3 about the existence and approximation of a minimal solution hold true in a 
general non Markovian framework with the following standard modifications : the terminal 
condition g^Xx) is replaced by a square integrable random variable ^ € L^(fi,.FT), the 
generator is a map / from x [0, T] x R x M*^ into R, satisfying a uniform Lipschitz 
condition in {y, z), and /(•, y, z) G L|(0, T) for all (y, z) € M x M*^, and the cost coefficient 
is a map c from Q x [0, T] x M x M*^ x into M, satisfying a uniform Lipschitz condition in 
[y, z), and c{-,y, z, e) G L|(0, T) for all {y, z,e) eRxR'^ x E. 
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Remark 2.2 Without the /i-constraint condition (j2.1ip on jumps, we have existence and 
uniqueness of a solution (Y, Z,U, K) with = to (j2.10p . from results on BSDE with 
jumps in [21] and [2]. Here, under (j2.1ip on jumps, it is not possible in general to have 
equality in (j2.10p with K = 0, and as usual in the BSDE literature with constraint, we 
consider a nondecreasing process K to have more freedom. The problem is then to find a 
minimal solution to this constrained BSDE, and the nondecreasing condition (|2.8p on h is 
crucial for stating comparison principles needed in the penalization approach. The primary 
example of constraint function is h{u, e) = —u, i.e. nonpositive jumps constraint, which is 
actually equivalent to consider minimal solution to BSDE (|1.13p as showed later. 

2.2 The case of nonpositive jump constraint 

Let us recall the BSDE defined in the introduction: find a triplet (Y, Z, K) G 5^ x L^(W) x 
such that 

Yt = g{XT) + f{Xs,Ys,Zs)ds + KT-Kt (2.12) 
{Zs,dWs)+ [ [ c{X,-,Y,-,Zs,e)fi{ds,de), 0<t<T,a.s. 

Jt J E 

such that for any other triplet (F, Z, ^) G 5^ x L^(W) x satisfying (|2.12p . it holds that 

Yt < Yt, 0<t<T, a.s. 

We will call such Y (and, by a slight abuse of notation, {Y, Z,K)) the minimal solution to 
()2.12p . We claim that this problem is actually equivalent to problem (j2.10p - (j2.1ip in the 
case h{u, e) = —u, corresponding to nonpositive jump constraint condition: 

Ut{e) < 0, dP (g> dt (E) X{de) a.e. (2.13) 

Indeed, let {Y, Z, U, K) be any solution of (j2.10p and (|2.13p . Define a process K by dKt = 
dKt — J^Us{e)fi{dt,de), < t < T, then K is nondecreasing, and the triplet {Y,Z,K) 
satisfies (j2.12p . It follows that the minimal solution to (j2.12p is smaller than the minimal 
solution to (j2.10p and (j2.13p . We shall see in the next section, by using comparison principles 
and penalization approach, that equality holds, i.e. 

minimal solution Y to ([2TT2]) = minimal solution Y to (pTOj) . (pT3]) . 




We shall illustrate this result by considering a special: when the functions / and c do 
not depend on y, z (i.e., the impulse control case). In this case, one can obtain directly the 
existence of a minimal solution to ()2.10p - ()2.13|) and (I2.12p by duality methods involving 
the following set of probability measures. Let V be the set of V <^ if-measurable essentially 
bounded processes valued in (0,oo), and given u £ V, consider the probability measure P'^ 
equivalent to P on {^},J^t) with Radon-Nikodym density : 



dP" 

If 



£t[J^ J^'^tie) - l)/2(dt,de)), (2.14) 
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where £t{-) is the Doleans-Dade exponential. Notice that the Brownian motion W remains 
a Brownian motion under P*^, which can then be interpreted as an equivalent martingale 
measure for the "asset" price process W . The effect of the probability measure P'', by 
Girsanov's theorem, is to change the compensator \{de)dt of /i under P to ut{€)^{de)dt 
under Y"^ . 

In order to ensure that the problem is well-defined, we need to assume : 

(HI) There exists a triple {Y , Z, K) e x L'^(W) x satisfying (|2T^ . 

This assumption is standard and natural in the literature on BSDE with constraints, and 
means equivalently here (when / and c do not depend on y, z) that one can find some 
constant y S ffi, and Z G L^(W) such that 

y+ riZs,dWs) > g{XT)+ r f{Xs)ds+ r [ c{X,-,e)fi{ds,de) a.s. 
Jo Jo Jo Je 

This equivalency can be proved by same arguments as in [8]. Notice that Assumption (HI) 
may be not satisfied as shown in Remark 13.11 in which case the problem (|2.12p is ill-posed. 

Theorem 2.1 Suppose that f and c do not depend on y,z, and (HI) holds. Then, there 
exists a unique minimal solution {Y, Z, K,U) G 5^ x L^(W) x L^(/2) x A^, with K pre- 
dictable, to (I2.10|) - (j2.13p . Moreover, {Y,Z,K) is the unique minimal solution to (|2.12p with 
Kt = Kt — Jq Us{e) fi{ds, de) , and Y has the explicit functional representation : 



Yt = ess sup W 



g{XT)+ r f{Xs)ds+ r I ciX,-,e)fiids,de) 
Jt Jt Je 



, 0<t<T. 



Proof. First, observe that for any (y,Z, [7,^) € S^xL^{W)xL'^{jl)xA'^ (resp. {Y,Z,K) 
€ X L2(W) X A^) satisfying ([230]) - (|2T3D (resp. (|2T2]) ). the process 

Qt ■■= Yt+ [ f{Xs)ds+ [ [ c{X,-,e)fi{ds,de), < t < T, 
Jo Jo Je 

is a P'^-supermartingale, for all G V, where the probability measure P'^ was defined in 
(fmp . Indeed, from (f230]) - (l233D (resp. (|2T2]) 1. we have 



Qt = Qo+ I {Zs,dWs)-Kt, with Kt = Kt- [ Us{e)fi{ds,de), 
Jo Jo 

( resp. Qt = Qo+ [ {Zs,dWs) - Kt), 0<t< T. 

Jo 

Now, by Girsanov's theorem, W remains a Brownian motion under P'^ , while from the 
boundedness of G V, the density dP'^ /dF lies in L^(P). Hence, from Cauchy-Schwarz 
inequality, the condition Z G L^(W), and Burkholder-Davis-Gundy inequality, we get the 
P'^-martingale property of the stochastic integral J {Z, dW), and so the P^'-supermartingale 
property of Q since K (resp. K) is nondecr easing. This implies 



Yt > 



Yt+ I f{Xs)ds+ f [ c{X,-,e)fi{ds,de) J't 
Jt Jt Je 
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and thereby, from the arbitrariness of P'^, v £ V, and since Yt = giXj-) 



Yt := ess sup W giXr) + / f{Xs)ds + 
i/GV l 



c{X,-,e)^l{ds,de) < Yt.{2.15) 



To show the converse, let us consider the process Y defined in (j2.15p . By standard 
arguments as in [11], the process Y can be considered in its cad-lag modification, and we 
also notice that Y £ S'^. Indeed, by observing that the choice of z/ = 1 corresponds to 
the probability P"" = P, we have Y <Y <Y, where {Y, Z,K) e x L^(W) x is a 
solution to (j2.12p . and 



Y 



E 



giXr) + / f{Xs)ds + 



T 



c{Xg- , e)fj,{ds, de) 



Thus, since Y lies in S"^ from the linear growth conditions on (7, /, and c, and the estimate 
(|2.2|) . we deduce that y G 5^. Now, by similar dynamic programming arguments as in [llj . 
we see that the process 



Qt = Yt+ f f{Xs)ds+ f [ c{X,-,e)fi{ds,de), 0<t<T, 
Jo Jo Je 



(2.16) 



lies in 5^, and is a P'^-supermartingale, for all z/ G V. Then, from the Doob- Meyer decom- 
position of Q under each P'', z/ S V, we obtain : 



Qt 



Yo + M" - K", 



(2.17) 



where M'^ is a P'^-martingale, = 0, and K'^ is a P'^ nondecreasing predictable cad-lag 
process with Kq = 0. Recalling that is a P'^-Brownian motion, and since fjy{ds,de) 
:= fi{ds,de) — Vs[e)\{de)ds is the compensated measure of /.i under P"^, the martingale 
representation theorem for each M'^ , u gives the existence of predictable processes Z'^ 
and such that 

Qt = Yo+ [ {Z^,dWs)+ [ [ U^{e) fluids, de)-K^, 0<t<T. (2.18) 
Jo Jo Je 

By comparing the decomposition (|2.18|) under P'^ and P corresponding to = 1, and 
identifying the martingale parts and the predictable finite variation parts, we obtain that 
= Z^ =: Z, U" = U^ =:U for all z^ G V, and 



= K} - f [ Us{e){us{e) - l)X{de)ds, 0<t<T. 
Jo Je 



(2.19) 



Now, by writing the relation (|2.18|) with v = e > 0, substituting the definition of Q in 
(|2.16p . and since Yt = g{XT), we obtain : 



Yt = giXT)+ r f{Xs)ds- r {Zs,dWs)- r / 

Jt Jt Jt Je 

+ [ [ Us{e)e\{de)ds + K^-Kl, 0<t<T. 
Jt Je 



(Usie) - c{Xg- , e))n{ds, de) 



(2.20) 



10 



From (j2.19|) . the process K"^ has a hmit as e goes to zero, which is equal to = 
+ Jq Us{e)X{de)ds, and inherits from K^, the nondecreasing path and predictabihty 
properties. Moreover, since Q G 5^, in the decomposition ()2.17|) of Q under P = P'^ for 
= 1, the process lies in and e A^. This imphes that Z e L2(W), [/ G L2(/2), 
and also that G A^. By sending e to zero into ^I7IU\f . we obtain that {Y,Z,U,K^) € 
X L^(W) X L^(/2) X A^ is a solution to (|2.1U|) . Let us finally check that U satisfies the 
constraint : 

Ut{e) < 0, dP (g)dt(g) X{de). (2.21) 

We argue by contradiction by assuming that the set F = {{uj,t,e) € ^lx[0,T]x E : Ut{e) > 
0} has a strictly positive measure for dP x dt x X{de). For any A; > 0, consider the process 
Vk = ^F'= + {k + 1)1f, which lies in V. From (j2.19p . we have 



E[K^'=] = B[K^]-kB\ [ [ lFUt{e)\{de)dt 



< 0, 



for k large enough. This contradicts the fact that K!^^ > 0, and so (j2.2ip is satisfied. 
Therefore {Y, Z,U, K^) is a solution to (|2.1U|) - (|2.13|) . and it is a minimal solution from 
(I2.15p . Y is unique by definition. The uniqueness of Z follows by identifying the Brownian 
parts and the finite variation parts, and the uniqueness of {U, K^) is obtained by identifying 
the predictable parts by recalling that the jumps of fi are inacessible. By denoting = 
— Jq Us{e)fi{ds, de), which lies in A'^, we see that {Y,Z,K°) is a solution to (l2T2]) . 
and it is minimal by (j2.15p . Uniqueness follows by identifying the Brownian parts and the 
finite variation parts. □ 

Remark 2.3 In Section HJ we shall relate rigorously the constrained BSDEs (I2.10j) - (l2.1ip 
to QVIs. In particular, the minimal solution Yt to (j2T0]) - (f2T3]) or (f2l^ is Yt = v{t,Xt) 
where v is the value function of the impulse control problem (jl.4p . Together with the 
functional representation of Y in Theorem 12. 1|, we then have the following relation at time 
t = 0: 

v{0,Xo) = supB''\g{XT) + r f{Xs)ds+ r [ c{X,-,e)ij{ds,de)\. (2.22) 
vgv l jo Jo Je 

We then recover a recent result obtained by Bouchard [1], who related impulse controls to 
stochastic target problems in the case of a finite set E. We may also interpret this result as 
follows. Recall that the effect of the probability measure P^ is to change the compensator 
X{de)dt of fj. under P to i't{e)X{de)dt under P'^. Hence, by taking the supremum over all 
P'^, we formally expect to retrieve in distribution law all the dynamics of the controlled 
process in (jl.Sp when varying the impulse controls a, which is confirmed by the equality 
(|2:22D . 



Finally, we mention that the above duality and martingale methods may be extended 
when the generator function / is linear in z by using Girsanov's transformation. Our main 
purpose is now to study the general case of /i-constraints on jumps, and nonlinear functions 
/ and c depending on y,z. 
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3 Existence and approximation by penalization 

In this section, we prove the existence of a minimal solution to (j2.10p - (|2.1ip . based on 
approximation via penalization. For each n G N, we introduce the penalized BSDE with 
jumps 

y," = g{XT)+ r f{Xs,Y^,Z^)ds + n r f h-{U^{e),e)\{de)ds (3.1) 

Jt Jt JE 

- r {Z'-,dWs)- r f {U^{e)-c{X,-,Y^^.,Z^,e))fi{ds,de), 0<t<T, 
Jt Jt Je 

where h~{u,e) = max(— /i(n, e), 0) is the negative part of the function h. Under the 
Lipschitz conditions on the coefficients /, c and h, we know from the theory of BS- 
DEs with jumps, see [21] and [2], that there exists a unique solution (y",Z",{7") S 
52 X L2(W) X L2(/2) to dsn). We define for each n G N, 

KJ" = n [ [ h-{U^{e),e)X{de)ds, 0<t<T, 
Jo Je 

which is a nondecreasing process in A"^. The rest of this section is devoted to the conver- 
gence of the sequence {Y"', Z", [/", K^)n to the minimal solution we are interested in. 

3.1 Comparison results 

We first state that the sequence (y")n is nondecreasing. This follows from a comparison 
theorem for BSDEs with jumps whose generator is of the form f{x,y,z,u) = f{x,y,z) + 
J^h{u{e),e)X{de) for some nondecreasing function h, which covers our situation from the 
nonincreasing condition on the constraint function h. 

Lemma 3.1 The sequence {Y"')n is nondecreasing, i.e. for all n E N, 1^" < Y^'^^ , < 
t < T, a.s. 

Proof. Define the sequence {y^)n oiV ® (^-measurable processes by 

Vr{e) = Urie) - ciXt- , Y^l , Zf, e), (t, e) G (0, T] x E and 
V^{e) = U^{e)-c{X,,Y^,Z^,e), eeE, 

From ()3.ip and recalling that X and Y are cad-lag, we see that Z", V^) is the unique 
solution in x L^(W) x L^(/i) of the BSDE with jumps : 

y," = g{XT)+ r Fn{Xs,Y:',Z2,Vr)ds- [^{Z^^dW^)- C f Vr{e)il{ds,de), 
Jt Jt Jt Je 

with En{x,y, z,v) = f{x,y,z) + J^(v{e) + nh~{v{e) + c{x,y, z,e),e))\{de). Since h~ is 
nondecreasing, we have 



Fn{t,x,y,z,v) - Fn{t,x,y,z,v') = / {{v{e) - v' (e)) + n[h {v{e) + c{x,y, z,e),e) 

Je 

- h~{v'{e) +c{x,y,z,e),e)]} \{de) 
- I {0- + '^{v{e)>v'{e)}nkh){v{e) -v'{e))] \{de). 



E 
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Moreover, since Fn+i > Fn, we can apply the comparison theorem 2.5 of [20], and obtain 
that y^" < yi''+\ < t < r, a.s. □ 

The next result shows that the sequence (Y^)n is upper-bounded by any solution to the 
constrained BSDE. Arguments in the proof involve suitable change of probability measures 
P", u £V, introduced in ([231]) . 

Lemma 3.2 For any quadruple {Y, Z, U, K) £ S"^ x L^(W) x L^ijl) x satisfying ([230]) - 
(j2.1ip . and for all n G N, we have 

Yt < 0<t<T, a.s. (3.2) 

Moreover, in the case : h{u, e) = —u, the inequality (j3.2p also holds for any triple (Y , K) 
G 5^ X L2(W) X A^ satisfying (|2T2]) . 

Proof. We state the proof for quadruple (Y^Z^U^K) satisfying (j2.10|) - (j2.1ip . Same argu- 
ments are used in the case : h{u,e) = —u and {Y,Z,K) G 5^ x L^(W) x A^ satisfying 

Denote Y = Y-Y"", Z = Z-Z''\ f = f{X, Y, Z)-f{X, F", Z") and c = c{X- ,Y~,Z, e)- 
c{X - , y" , Z", e). Fix some v £ V (to be chosen later). We then have : 



Yt= [ fsds+ [ [ csti{ds,de)- [ {Z„dWs)- [ [ \ Us {e) - Ul'{e)} fluids, de) 
-/ / [Usie)-U2{e)\i^sie)X{de)ds-n f j h' {U2{e),e)X{de)ds + Kt - Kt, 



where fl'^{dt,de) = n{dt,de) — ut{e)\{de)dt denotes the compensated measure of /i under 
Y*^ . Let us then define the following adapted processes: 

f{XuYt,Zt)- f{Xt,Y(^,Zt) ^ 
- ^ ^{?i7^0}> 

and h the M^- valued process defined by its i-ih. components, i = 1, . . . , d: 

.i _ fiXt,Y,^,zt'^) - f{Xt,Yr,zP) ^ 

- yi -^{V^VO}' 

where Z^"^ is the M'^-valued random vector whose i first components are those of Z and 
whose {d — i) lasts are those of Z", and is the i-th component of z[^ — zf\ Let us 
also define the V ^ (^-measurable processes (5 in R and ^ in M'^ by: 

, c(x,-,y,-,zo-c(x,-,y-,Zi,e) 

1 J. 



and 



_ c(x,- , y,"„ , zr^' , e) - c(x,- , y,"_ , zr> , e) ^ 
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Notice that the processes a, b, 6 and £ are bounded by the Lipschitz conditions on / and c. 
Define also = at + j^5t{e)vt{e)X{de), (3^ = bt + ^((e)f( (e)A(de), which are bounded 
processes since a,b,6,i are bounded and A is a finite measure on E, and denote V"{e) = 
Ut{e) — Up{e) — 6t{e)Yt — £t{e) ■ Zt. With these notations, and recalhng that h~{Us{e)) = 
from the constraint condition (|2.11|) . we rewrite the BSDE for Y as: 



Yt 



rT 




t JE 



V^ie)fi''{ds, de) + Kt - Kt 



+ 



n[h-{Us{e),e) - h-{U^{e),e)] - Vs{e)[Us{e) - U^{e)]] \{de)ds. 



Consider now the positive process solution to the s.d.e.: 



dT'; 



Vdatdt + iP^t^dWt)), = 1, 



and notice that lies in 5^ from the boundeness condition on a'^ and j3'^ . By Ito's formula, 
we have 



dr^yt = [n[h-{iJt{e),e) - h-{Ur{e),e)] - Me)[Ut{e) - C/f (e)]} X{de)ds 

-V^dKt + T'iiZt, dWt) + T^tYt-iPu dWt) - / Vr{e)fi'{dt, de), 
which shows that the process 



n[h~{Us{e),e) - h-{U^{e),e)] - Me)[Us{e) - U^{e)]] \{de)ds 



is a P'^-supermartingale as soon as {Ut{e),e) — h (C/"(e), e)] — z^t(e)[[7t(e) — ^7"(e)] >0 
for all {t, e) G [0, T] x E, and so 



^V", j^[n[h-iUsie),e) - /i-(C/r(e), e)] - i^!{e)[Us{e) - U^{e)]] X{de)di 
Now, from the Lipschitz condition on h, we see that the process ly^ defined by 



n^-(^4e),e)-.-g.r(e),e)1 ^ ^^^^^ h-{U^{e), c) > 

e else 



is bounded and so lies in V, and therefore by taking = we obtain : 

rT r 
IE 

From Bayes formula, this is written as : 



rfYt > -eW 



s / [Us{e) ^r'(e)]l{[/^(e)>f/n(e)}u{h-(C/,"(e),e)=0}'^('^^)'^'^ 



r'^Yt 



> -eE 



(3.3) 



Zjrp 

Jt 



s / [t^«(e) ^r(e)]l{[7^(e)>;7n(g-)}u|;j-(f/n(e)^e)=o}'^('^^)'^'^ 



14 



where Z^^ is the Doleans-Dade exponential : 

= exp 
By definition of v'^ ^ we have 



y (z^f (e) — l)A(de)ds^ exp log{ul{e))fi{ds,de) 



< ^ exp yj^ {nkh - l)X{de)dsj exp 



{e - l)X{de)ds 



E 



where is the solution to dZ" = Zl'_ {nkh ~ 1) de), Zq = 1. This shows that 

2 

zf 



E 



is uniformly bounded for e in the neighborhood of 0"*^. Notice also that the family (r''^)o<e<n 
is uniformly bounded in 5^ so that 

2 



E 



/ [f^s(e) ^r(6)]l{(7s(e)>f/y(e)}U{/i-(C/,"(e),e)=0}-^(^'^)'^^ 



is again uniformly bounded for e in the neighborhood of 0"*". From the conditional Cauchy- 
Schwartz inequality, we deduce that 



E 



Zjr-p 



t Jt 



is uniformly bounded for e in the neighborhood of 0"*". Finally, since lim^^or^ = F^ > 
0, by sending e to zero into ()3.3p . we conclude that > 0. □ 



3.2 Convergence of the penalized BSDEs 

We impose the following analogue of Assumption (HI). 

(H2) There exists a quadruple {Y,Z,k,U) £ x L^{W) x L^{jl) x satisfying 
(|230D - (f2m) . 

Assumption (H2) ensures that the problem (j2.10p - (j2.1ip is well-posed. As indicated in 
paragraph 12.21 Assumption (H2) in the case h{u, e) = — u, is stronger than Assumption 
(HI). We provide in Section [B] some discussion and sufficient conditions under which (H2) 
holds. 

Remark 3.1 The following example shows that conditions (HI) and (H2) may be not 
satisfied : consider the BSDEs 



Y, 



I {Zs,dWs)+ I I c^i{ds,de) + KT-Kt, (3.4) 
Jt Jt Je 



and 



Yt = - J^iZs, dWs) - lEpsie) - c]fi{ds, de) + Kt - Kt 
-Us{e) > 



(3.5) 
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where c is a strictly positive constant, c > 0. Then, there does not exist any solution to 
(j3.4p or (j3.5p with component Y £ S"^. On the contrary, we would have 

Yo > - [ {Zs,dWs)+cfii[0,T]xE), a.s. 
Jo 

which implies that for all n G N*, = n G V, 

Yo > E'^ - / {Zs,dWs) +cn{[0,T] X E) = cnX{E)T. 
. Jo 

By sending n to infinity, we get the contradiction : ||1^||52 = oo. 

We now establish a priori estimates, uniform on n, on the sequence {Y", Z", U"^,K'^)n- 



Lemma 3.3 Under (H2) (or (HI) in the case : h{u,e) 
C such that 



-u), there exists some constant 



lY'^W 2 + ||Z"|| , + ||C/"|| , + \\K 

' II c2 I II III 2fwl I II IIT.2/-N I II 



■ II- ii,2,w) ' n,2„) I n^^ n,2 < C, Vn G N. 



(3.6) 



Proof. In what follows we shall denote C > to be a generic constant depending only on 
T, the coefficients /, c, the process X, and the bound for Y in (HI) or (H2), and which 
may vary from line to line. 

Applying Ito's formula to ji^'^P, and observing that K"^ is continuous and AY^ = 

U^{e) - c{X,- , yi , Z^, e)fi{{t}, de), we have 

E|5(Xr)P = E|y,"|2-2E C Y^ f{X,,Y^ , Z^)ds - 2^ C Y^dK^ + ^ ^ iZ^fds 

Jt Jt Jt 

+ B r [ {\Y:^^+U^{e)-ciX,^,Y:i,Z^,e)\'-\Y:i\'}X{de)ds 

Jt JE 

Prom the linear growth condition on / and the inequality Y^ < by Lemma 13.21 under 
(H2) (and also under (HI) in the case h{u,e) = —u), and using the inequality 2ab < 
+ ab^ for any constant a > 0, we have: 



Eiy.'^P + E 



|Z"P(is + E 



T 



|C/;(e) - c(X,- , y;_ , Z^, e)\'Xide)ds 
< B\g{XT)\^ + 2CB j^' \YP\ (1 + + |y,"| + \Z^\) ds 



-2E 



T 



1. 



Y^.iKie) - c{X,- , Y;^ , , e))X{de)ds + -E 



sup \Yt\ 



« ^te[o,T] 



n\2 
t I ■ 



Using again the inequality 2ab < j^a + ab , in particular for a = 2, yields 



E|y:"r + E 



|Z"|2(is + -E 
2 



T 



\U^{e) - c{X,-,Y^^^,Z^,e)\^X{de)ds 
< E|5(Xt)P +2CE^' lYJ'lil + lXsl + lY^Pl + lZ^Dds 



+ 2A(£;)E / |y"P(is + -E 

a 



< C(l + E / m'^ds] +Ib 



sup \Yt 

te[o,T] 

T 



n|2j„ I „.T?\T^n T^n\2 



\Z^\'ds + aE\K!^- 
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Then, by using the inequahty (a — b)^ > — b^, we get 



T 



\Z^\^ds + -E 



T 



\U^{e)\'^\{de)ds 



1, 



< -EJ^ jJc{X,-,Y:.,Z^,e)\'\{de)ds 
rT 

+ C(1 + E/ lY^'l^ds] +aE\K^ - Kr^l^ 



T 



< C{l + E I \Y:fds]+-El \Z^\'^ds + aE\K^- K^\'^, 



il 
4 

from the hnear growth condition on c. Now, from the relation 

rT 



(3.7) 



+ r [ {U:ie) - c{X,-,Yp.,Z2))fi{ds,de) + Z^.dW^, 

Jt JE Jt 

and the hnear growth condition on /, c, there exists some positive constant Ci s.t. 



n|2 
t I 



E\K^ - K 



< Ci{l + E\Yt"'\^ + E lYj'l^ + \Z^\'ds + E 



T rT 
"1 2 I TTi / w^n\2 I I ryn\2 



\U^{e)\'^\{de)ds] . (3.8) 
It Jt je J 

Hence, by choosing a > s.t. Cia < 1/4, and plugging into (j3.7|) . we get 

ElY^^l^ + E r \Z^\'^ds + E r [ \U^{e)\'^\{de)ds < C ( 1 + E T \Yp\^ds 
Jt Jt Je \ Jt 

By applying Gronwall's lemma to t ^ E|l^"p and p.Sp . we obtain 

sup E|yj"|2 + E / \Z2\'^ds + E [ [ \U^{e)\'^X{de)ds + E\K^\^ < C. (3.9) 
o<t<r JO Je 

Finally, by writing from ()3.ip that 



sup I y," I < \g{XT)\+ r \f{Xs,Ys,Z,)\ds + KJ^ 

<t<T Jo 

r [ \U2ie)-ciX,-,Y,-,Z,,eMds,de), 
Jo Je 









+ sup 


r {Zs,dWs) 




se[o,T] 


Jo 


Jo Je 



we obtain the required result from the Burkholder-Davis-Gundy inequality, the linear 
growth condition on /, c, and (|3.9p . □ 



Remark 3.2 A closer look at the proof leading to the estimate in (j3.6p shows that there ex- 
ists a universal constant C, depending only on T, and the linear growth condition constants 
of /, c, such that for each n £ N : 



sup E[yj"]2 < c(i + e|^(Xt)P + e 

tG[0,T] ^ 



r /"^ 1 








/ \Xt\^dt 


+ E 


sup \Yt\^ 


) 






^te[o,T] J 





(3.10) 
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Lemma 3.4 Under (H2) (or (HI) in the case : h{u,e) = —u), the sequence of processes 
iXl^) converges increasingly to a process (It) with Y ^ S"^ . The convergence also holds in 
L|(0, T) and for every stopping time r G [0, T], the sequence of random variables (Y^) 
converges to Yr in L^(r2,^T-), i.e. 



lim E 

n— »oo 



y^" - Yti'dt 



and lim E 

n— ►oo 



0. 



Proof. From Lemmas 13.11 and 13. 2| the (nondecreasing) limit 



Y, 



lim y", < t < T, 



(3.11) 



(3.12) 



exists almost surely, and this defines an adapted process Y. Moreover, by Lemma 13.31 and 
convergence monotone theorem, we have 



E 



sup \Yt\ 

0<t<T 



< oo. 



From the dominated convergence theorem, we also get the convergences ()3.1ip . It remains 
to check that the process Y has a cadlag modification. We first show that (y")n are quasi- 
martingales with uniformly bounded conditional variations. That is, there exists a constant 
C such that, for any partition vr : = to < < " " " < = 

n-l 

E{|yp| + ^ mY,l^\TtA - Yul} < C, Vtt, Vn. (3.13) 

1=0 

In fact, by (|3.ip we have 

n— 1 n— 1 ft +1 

e{ ^ \E{Y,l^\J^tA - Y,^\} = e{ ^ |e[ / f{Xs, Yr, Z^)ds 

+n r""' / h-{U:{e),e)X{de)ds- T^' [ {U^ie) - c{X,- ,Yr. , ,e))X{de)ds J^J } 
Jti Je Jti Je ^ ^ 



< E 



\fiX,,Yr,Z^)\ds + 







\U2ie) - c{X,-X:-.Z^,e)\\{de)ds + KJ]^ 



'0 JE 

Recah (I23D, ([23]), and ([22]), we have 

n-l 



^[\Y^\+Y,\nyuj^u} -yu\] 

i=0 

< Ce{i + \Xt\ + j\l + \Xs\ + + \Z':\]ds + £ j \K{e)\\{de)ds + ^C?}. 

Applying (j2.2p and Lemma 13.31 we obtain (|3.13p immediately. Now by Meyer-Zheng [15] 
(or see |14]). there exists a subsequence (y"'=)fc and a cadlag process y such that (y"'=)fc 
converges to Y in distribution. On the other hand, by (j3.12p . (y"''')fc converges to Y, P-a.s. 
Then Y and Y have the same distribution, and thus Y is also cadlag. □ 

We now focus on the convergence of the difi'usion and jump components {Z"^, U'""). In 
our context, we cannot prove the strong convergence of {Z"^, C/") in L^(W) x L^(/2), and so 
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the strong convergence of j^Z'^dW and /o /^C/"(s,e)/i(ds,de) in L^{n,J^t), see Remark 
13. 3i Instead, we follow and extend arguments of Peng [18], and we shall prove that (Z", [/") 
converge in LP(W) x LP(/i), for 1 < p < 2. First, we show the following weak convergence 
and decomposition result. 

Lemma 3.5 Under (H2) (or (HI) in the case: h{u,e) = —u), there exist (j) G L|(0, T), 
Z G L^(W), V G L^(/i) and K G predictable, such that the limit Y in \?>.\2\ has the 

form, 

Yt = Yo- f ct)sds-Kt+ [\zs,dWs)+ [ [ Vs{e)fi{ds,de), 0<t<T.{3.U) 

Jo Jo Jo J E 

Moreover, in the above decomposition of Y , the components Z and V are unique, and 
are respectively the weak limits of (Z^) in L^(//) and of {V") in L^(/i) where Vt"'{e) = 
C/"(e) — c(Xt- ,Y_^, Z^,e), (j) is the weak limit in L^(0,T) of a subsequence of (/") := 
(/(X, y", Z^)), and K is the weak limit in Lp(0, T) of a subsequence of {K^)- 



Proof. By Lemma 13.31 a^^id the linear growth conditions on /, c together with (j2.2p . 
the sequences (/"), {Z'^), {V^) are weakly compact, respectively in L^(0,T), L^(W) and 
L^(/2). Then, up to a subsequence, (/"), {Z^), {V^) converge weakly to Z and V. By 
Ito representation of martingales, we then get the following weak convergence in L^(ri,.7v) 
for each stopping time r < T: 



r f^ds - r (t>sds, r {z^,dWs) ^ r{zs,dws 

Jo Jo Jo Jo 



Vj'{e)n{ds,de) / / Vs{e) ^x{ds , de) . 

E Jo Je 



Since, we have from p.ip : 



= -y;^ + yo"- r fys+ r{Z^,dWs)+ [ [ Vr{e)f,{ds,de), (3.15) 
Jo Jo Jo Je 

we also have the weak convergence in L^(r2,.7v) : 

^ Kr := -Yr + Yo - [ 4>sds+ I {Z„dWs)+ [ [ y,(e)/i((is, de). (3.16) 

Jo Jo Jo Je 

The process K inherits from K"' the nondecreasing path property, is square integrable, 
cad-lag and adapted from (j3.16p . and so lies in A^. Moreover, by dominated convergence 
theorem, we see that converges weakly to K in L'^(0,T). Since is continuous, and 
so predictable, we deduce that K is also predictable, and we obtain the decomposition 
()3.14p for Y. The uniqueness of Z follows by identifying the Brownian parts and finite 
variation parts, and the uniqueness of V is then obtained by identifying the predictable 
parts and by recalling that the jumps of /U are inacessible. We conclude that (Z, V) is 
uniquely determined in (j3.14p . and thus the whole sequence (Z",y") converges weakly to 
(Z,y) in L2(W) X L2(/2). □ 

The sequence (f/") is bounded in L^(/2), and so, up to a subsequence, converges weakly 
to some U G L^(/2). The next step is to show that the whole sequence (C/") converges 
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to U and to identify in the decomposition (j3.14p cpt with f{Xt,Yt,Zt), and Vt{e) with 
Ut{e) — c{X^- , Y^- ,Zt,e). Since / and c are nonUnear, we need a resuh of strong convergence 
for (Z^) and (f/") to enable us to pass the hmit in f(Xt,Yf^,Zf}) as weh as in C/"(e) — 
c(Xj- , Y^ , Z", e), and to eventually prove the convergence of the penalized BSDEs to the 
minimal solution of our jump-constrained BSDE. We shall borrow a useful technique of 
Peng [IH] to carry out this task. 

Theorem 3.1 Under (H2), there exists a unique minimal solution (Y, Z,U, K) G 5^ x 
L2(W) X L2(/2) X with K predictable, to (f230|) - (pTT]) . Y is the increasing limit of (Y'^) 
in (j3.12p and also in L|(0, T), K is the weak limit of (K^) in L|(0, T), and for any p G 
[1,2), 

11-^" - ^IIlp(w) + 11^" ~ ^IIlp(m) * 0' 

as n goes to infinity. Moreover, in the case : h{u, e) = —u, (y, Z, K) is the unique minimal 
solution to (j2.12p with Kt = Kt — f^Us{e)^{ds,de), and this holds true under (HI). 
Consequently, the minimal solution Y to (|2.12p and to (|2.10p - ()2.13p are the same. 

Proof. We apply Ito's formula to jl^" — Itp on a subinterval {a, r], with < a < t < T, 
two stopping times. Recall the decomposition (|3.14p . (|3.15p of Y, y", and observe that if" 
is continuous, and A(yi" - Yt) = AKt + /g;(V;"(e) - Vs{e))fi{{t},de). We then have : 

E|y;^ - y,|2 = e\y^' - y^\^ + e T \z^ - Zs\^ds + 2e T^" - y][(A. - f^ds 

- 2E /"[y," - y,]dJC + 2E / [y^"L - Y,-\dK, + E ^ \AKtf 

J <J J (cr,Tl 1 

+ E / / [|y« _y^_+y-(e)-V;(e)|2-|l^'i -y,-|V(rfs,de) 

J{o-,t] 

= E|yj^ - y,|2 + E r - z.fds + 2E - Y,\{4>, - f^]ds 

-2E I K" - Ys]dK^ + 2E / [y^"_ - y,- + AKs]dKs 

J (J J (cr,r] 

-E ^ \AKt\^ + Er I \Vr{e)-Vs{e)\^X{de)ds 

+ 2E r / (y,"-y,)(y,"(e)-y,(e))A(de)ds. 

Since (y^"-y,)dK^" < 0, and by using the inequality 2ab > -^-2b'^ with a = VJ'{e)-Vs{e) 
and b = Y^ — Yg, we obtain : 

E 1^ \Z^ - Z,\^ds + ^^1^1^ \Vr{e) - Vs{e)fds 

< B\Y^ - Yr\^ + 2E r |y," - Ysl^ds + 2E T |y," - y,||0, - f^\ds 

J a J cr 

i / \Y^--Y,- + AKs\dK, + EY,\^Kt\^- (3-17) 



2E 
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The two first terms of the right side of (j3.17p converge to zero by (j3.11|) in Lemma [3.4[ The 
third term also tends to zero since {(p — f^)n is bounded in L^(0,T), and so by Cauchy- 
Schwarz inequahty: 

B£\Y,^-Ys\\<t>s-mds < C(B£\YP-Ys\^dsy ^ 0. (3.18) 

For the fourth term, we notice that the jumps of y" are inacessible since they are determined 
by the Poisson random measure /i. Thus, the predictable projection of Y^ is ^1^" = Y^ . 
Similarly, from (j3.14p . and since K is predictable, we see that ^Yj = Y^- — AKt- Since 
yn increasingly converges to Y, then ^Y^ also increasingly converges to ^Y, and by the 
dominated convergence theorem, we obtain: 



lim E / \Y^^ -Y,- + AKs\dKs = 0. (3.19) 

'{0,T] 



For the last term in (|3.17|) . we apply Lemma 2.3 in [TS] to the predictable nondecreasing 
process K: for any 6,e > 0, there exist a finite number of pairs of stopping times {(Jk,Tk), 
k = 0, . . . ,N, with < (jfc < Tfc < r, such that all the intervals (cr^, r^] are disjoint and 

N N r- 

EY,{r, - a,) > T - eJ] J] {AKt)'<j. (3.20) 

fc=0 fc=Oo-fe<t<rfc 



We should note that in [18] the filtration is Brownian, therefore it is continuous, and hence 
each stopping time 0"^ can be approximated by a sequence of announceable stopping times 
In our case the stopping times a^s are constructed as the successive times of jumps of the 
predictable process K with size bigger than some given positive level, the approximation of 
(Tfc by announceable stopping times is again possible. We can thus argue exactly the same 
way as in Lemma 2.3 in |18j to derive both estimates in (j3.20p . 

We now apply estimate (|3.17p for each a = ak and r = r^, and then take the sum over 
k = 0,...,N. It follows that 



5;E r\Z--Zfds + lY.^ r [ \Vr{e)-Vs{e)\'ds 
< y2E\Y,l-Yrf + 2E |y--y,|2ds + 2E / \Yp - Y,\\<Ps - f^ds 

k=Q ^0 Jo 

N 

+ 2E / |y;i - Y,- + AKs\dKs + V E V \AKt\'^ . 

J(0,T] , 



JV 

\Y^''.-Y,- + AKs\dKs + 

'(O'^l k=0 t&(ak,Tk, 

From the convergence results in Lemma 13.41 ()3.18p and (j3.19p . we deduce that 



limsupVE /" ' - Z,|2ds + ^ Ve f " I iVj'ie) -Vs{e)\'^X{de)ds 



2 < £^ 

3 ' 



fc=0 tG{(7fc,rfc] 
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Thus, there exists an integer > such that for all n > l^s, we have 



A,'=0 

This implies 
and 



k=0 



o-fc J E 



V^{e)-Vs{e)\^\{de)ds < 



N 



(s,cj) e |J(cTfcH,rfcH] X n : \Z2{oj) - Zs{lo)\^ > 6 



k=0 



N 



< 



{s,e,Lo) G \J{ak{uj),n{uj)] xQxE : \Vj'{e,Lo) - Vs{e,oo)\^ > 5 



k=0 



< e. 



Together with ([320]), it follows that dt(g)P [(s, cu) G [0, T] x : - Zs{uj)\^ > 5] < e, 

and 

dt®XxP[{s,e,uj)£[0,T]xExn ■.\VP{e,uj)-Vs{e,uj)\'^>6] < e{l + X{E)). 
We deduce that for all 6 > 0, 

lim dt(g)P\{s,uj) £[0,T]xn -.iZ^iu)- Z,{u)\'^ >5] = 

and 

lim dt<S) X<S)P\{s,e,uj) £ [0,T]x E X n ■.\Vj'{e,uj) -Vsie,uj)\'^ > 5] = 0. 

This means that the sequences {Z^)n and (y'^)n converge in measure respectively to Z 
and V. Since they are bounded respectively in L'^(W) and L'^(/2), they are uniformly 
integrable in LP(W) and LP(/i) for any p G [1,2), respectively. Thus, (^") and (V*) 
converge strongly to Z and V in LP(W) and LP(/2), respectively. Recalling that J7"(e) 
= Vj"(e) + c(Xj- , y^" , , e), and by the Lipschitz condition on c, we deduce that the 
sequence (U^) converges strongly in LP(/i), for p G [1,2), to U defined by : 



Ut{e) 



Vtie) + c{Xt- , Yt- ,Zt,e), 0<t<T, e£ E. 



By the Lipschitz condition on /, we also have the strong convergence in Lp(0,T) of (/") 
= (/(X,y",Z")) to fiX,Y,Z). Since (/) is the weak limit of (/") in L|(0,T) , we deduce 
that (/) = f{X,Y,Z). Therefore, with the decomposition ()3.14|) and since Yr = lim„y^ = 
g{XT), we obtain immediately that (Y, Z, U, K) satisfies the BSDE (j2.10p . Moreover, from 
the strong convergence in L-'-(/I) of (C/") to U, and the Lipschitz condition on h, we have 



E 



/i~([/;(e),e)A(de)ds 



E 



h {Us{e),e)\{de)ds, 



as n goes to infinity. Since Kl} = n/g Jj^h-{U^{e),e)X{de)ds is bounded in l,^{fl,J^T) 
this implies 

fT 



E 



h {Us{e),e)X{de)ds 



0, 
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and so the constraint (j2.11|) is satisfied. Hence, (Y, Z, K, U) is a solution to the constrained 
BSDE (j2.10p - (j2.1ip . and by Lemma 13.21 Y = limY" is the minimal solution. The unique- 
ness of Z follows by identifying the Brownian parts and the finite variation parts, and then 
the uniqueness of {U, K) is obtained by identifying the predictable parts and by recalling 
that the jumps of jj, are inacessible. 

Finally, in the case h(u, e) = — n, the process 

Kt = Kt- [ [ Us{e)fi{ds,de), 0<t<T, 
Jo Je 

lies in A'^, and the triple (Y, Z,K) is solution to ()2.12p . Again, by Lemma 13.21 this shows 
that Y is the minimal solution to (j2.10p and to (j2.12p . The uniqueness of {Y,Z,K) is 
immediate by identifying the Brownian part and the finite variation part. □ 



Remark 3.3 From the estimate ()3.17p . it is clear that once the process K is continuous, 
i.e. AKt = 0, then {Z^,U"^) converges strongly to {Z,U) in L^(W) x L^(/i). This occurs 
in reflected BSDE's as in [10] or [12], see also Remark 14.31 In the case of constraints on 
jump component U as in (|2.1Up - (|2.1ip . the situation is more complicated, and the process 
K is in general only predictable. The same feature also occurs for constraints on Z as in 
|18j . To overcome this difficulty, we use the estimations (j3.20p of the contribution of the 
jumps of K, which allow to obtain the strong convergence of {Z^, U'"') in LP(W) x LP(/i) 
for p £ [1,2). Finally, notice that for the minimal solution (Y,Z,K) to the BSDE ([2T2]) . 
the process K is not predictable. 



3.3 The case of impulse control 



In the impulse control case (i.e. / and c depend only on X and h{u,e) = —u), we have 
seen in Theorem 12.11 that the minimal solution to our constrained BSDE has the following 
functional explicit representation : 



Yt = esssupE'^ 5(Xt)+ / f{Xs)ds + 



T 



c{Xg~ , e)fi{ds, de) 



In this case, we also have a functional explicit representation of the solution Y'"" to the 
penalized BSDE : 



Yr 



ess sup E'' 



T 



giXr) + / f{Xs)ds + 



c{Xg- )n{ds, de) 



(3.21) 



where V„ = {z^ G V ; z^s(e) < n V(s,e) € [0, T] x E a.s.}. Indeed, denote by y" the right 
side of ([Lll). By writing that (y",2"",f/") is the solution of the penalized BSDE (pTT]) . 
taking the expectation under P'^, for v G Vn, and recalling that is a P'^-Brownian 
motion, and v\{de) is the intensity measure of [i under P*^, we obtain : 



Yr 



giXr) + / f{Xs)ds + 



T 



c{Xs- , e)^i{ds, de) 



{n[C/;(e)]+-z.,(e)C/:(e)}A(de)ds 



(3.22) 
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Since this equality holds for any v G Vn, and observing that n[i7"(e)]+ — i/s(e)C/"(e) > 0, 
for all I' E Vn, we have 



T 



{n[;7,"(e)]+ - i/,(e)C/,"(e)} A((ie)ds 



Let us now consider the family {y^)e of V„ defined by 

' vmm± if [/n(e) > 
otherwise. 



(3.23) 



e 



Then, by using the same argument as in the proof of Lemma 13.21 we show that 



{n[^7,"(e)]+ - i/,(e)C/,"(e)} \(de)ds 



as e ^ 0, 



which proves with ^JIS^ that = . 

The representation (|3.21|) has a nice interpretation. It means that the value function of 
an impulse control problem can be approximated by the value function of the same impulse 
control problem but with strategies whose numbers of orders are bounded on average by 
nT\(E). This has to be compared with the classical approximation by iterated optimal 
stopping problems, where the n-th iteration corresponds to the value of the same impulse 
control problem but where the number of orders is smaller than n. The numerical advantage 
of the penalized approximation is that it does not require iterations. 



4 Relation with quasi-variational inequalities 

In this section, we show that minimal solutions to the jump-constrained BSDEs provide a 
probabilistic representation of solutions to parabolic QVIs of the form: 



mm 



dv 

— - Cv - f{-,v, a''D^v) , inf h{n''v - v, e) 

Ot e&E 



0, on [0,r) X M'^, (4.1) 



where C is the second order local operator 



1 



£v{t,x) = {b{x), Dxv{t,x)) — -tr{aa'' {x)D'^v{t, x)), 

and TC^, e £ E, are the nonlocal operators 

7i'^v{t, x) = v{t, X + 7(3;, e)) + c{x, v{t, x), a'' {x)Dxv(t, x), e). 

For such nonlocal operators, we denote for g G M'^ : 

H^[t,x,q,v] = v{t,x + ^{x,e)) + c{x,v{t,x),a'' {x)q,e). 

Note that when h{u) does not depend on e, and since it is nonincreasing in u, the QVI 
()4.ip may be written equivalently in 



mm 



dv 

— - Cv - f{-,v, a''D^v) , h{nv - v) 



0, on [0, T) X 
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with Hv = supgg^; W^t;. In particular, this includes the case of QVI associated to impulse 
controls for h{u) = —u, and /, c independent of y, z. 

We shall use the penalized parabolic integral partial differential equation (IPDE) asso- 
ciated to the penalized BSDE ()3.ip . for each n G N: 

_ ^ 
dt 



Cvn - f{-,Vn,cr''DxVn) - n h {n^Vn - Vn,e)X{de) = 0, (4.2) 

Je 



on [0, T) X M'^. 

To complete the PDE characterization of the function v, we need to provide a suitable 
boundary condition. In general, we can not expect to have v{T~ , .) = g, and we shall 
consider the relaxed boundary condition given by the equation: 



mm 



viT~, ■)- g , inf h(n^v(T~,-) -v(T-,-),e) = 0, on R'^, (4.3) 

e£E 



In the sequel, we shall assume in addition to the conditions of paragraph 12.11 that the 
functions 7, /, c, and h are continuous with respect to all their arguments. 



4.1 Viscosity properties 

Solutions of (j4.ip . (j4.2p and (j4.3p are considered in the (discontinuous) viscosity sense, and 
it will be convenient in the sequel to define the notion of viscosity solutions in terms of 
sub- and super-jets. For a locally bounded function u on [0, T] X we define its lower 
semicontinuous (Isc in short) u^,, and upper semicontinuous (use in short) envelope u* by 

u^{t,x) = liminf u(t',x'), u*{t,x) = limsup u(t',x'). 

{t',x')^{t,x),t'<T {t',x')->{t,x),t'<T 

Definition 4.1 (Suhjets and superjets) 

(i) For a function u : [0,T] x M'^ — > M, Isc (resp. use), we denote by J~u{t,x) the parabolic 
subjet (resp. J~^u{t,x) the parabolic superjet) ofu at {t,x) G [0,r] xM"^, as the set of triples 
{p, q,M) £RxR'^ xS'^ satisfying 

u{t',x') > {resp. <) u{t,x) +p{t' -t) + {q,x' - x) + ^{x' - x,M{x' - x)) 

+ o(|t' - t| + \x' - x|^). 

(a) For a function u : [0, T) xM.'^ ^ M., Isc (resp. use), we denote by J~u{t, x) the parabolic 
limiting subjet (resp. J~^u{t,x) the parabolic limiting superjet) ofu at {t,x) G [0,T] x M"^, 
as the set of triples {p, q, M) G M x M*^ x S'^ such that 

{p,q,M) = lim{pn,qn,Mn), {t,x) = lim(t„,x„), 

n n 

with {pn,qn,Mn) £ J~u{tn,Xn) (resp. J~^u{tn,Xn)), u{t,x) = limn(t„,x„). 

n 

We now give the definition of viscosity solutions to (j4.ip . (j4.2p and (j4.3p . 
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Definition 4.2 (Viscosity solutions to (j4.ip ) 

(i) A function u, Isc (resp. use) on [0,T) X W^, is called a viscosity supersolution (resp. 
subsolution) to (j4.ip if for each {t,x) e [0,r) x R'^, and any {p,q,M) G J~u{t,x) (resp. 
J~^u(t,x)), we have 



mm 



P- {Kx),q) - ^tT{aa''{x)M) - f{x,u{t,x),a''{x)q) 



inf h{7i'^[t,x,q,u] — u{t,x),e) > ( resp. <) 0. 

eei? J 

(ii) A locally bounded function on [0, T) x is called a viscosity solution to (|4.ip if is 
a viscosity supersolution and u* is a viscosity subsolution to (|4.ip . 

Definition 4.3 (Viscosity solutions to (j4.2p ) 

(i) A function u, Isc (resp. use) on [0,T) x M*^, is called a viscosity supersolution (resp. 
subsolution) to (|4.2p if for each {t,x) G [0,T) x M.'^, and any {p,q,M) E J~u{t,x) (resp. 
J'^u{t,x)), we have 

-p-{b{x),q) - ^tT{aa'{x)M) - f{x,u{t,x),a'{x)q) 

—n 1 h~{7i^[t,x,q,u] — u{t,x),e)X{de) > { resp. <) 0. 
JE 

(a) A locally bounded function u on [0,T) X is called a viscosity solution to (j4.2p if 
is a viscosity supersolution and u* is a viscosity subsolution to (|4.2p . 



Definition 4.4 (Viscosity solutions to (j4.3p ) 

(i) A function u, Isc (resp. use) on [0,T] x M*^, is called a viscosity supersolution (resp. 
subsolution) to (|4.3p if for each {x G W^, and any {p,q,M) G J~u{T,x) (resp. J'^u{T,x)), 
we have 



mm 



u{T, x) — g{x), inf h{7i^[T, x, q, u] — u(T, x), e) 



> ( resp. <) 0. 



(ii) A locally bounded function u on [0, T] x R'^ is called a viscosity solution to ()4.3p if 
is a viscosity supersolution and u* is a viscosity subsolution to (|4.3p . 



Remark 4.1 An equivalent definition of viscosity super and subsolution to (|4.3p . which 
shall be used later, is the following in terms of test functions : a function u, Isc (resp. use) 
on [0,r]x]R^is called a viscosity supersolution (resp. subsolution) to (14. 3p if for each (t, x) 
S [0, T) X W^, and any ip € C^'^([0, T] x M.'^) such that {t, x) is a minimum (resp. maximum) 
global oi u — ip, we have 



mm 



u(T,x) — g{x), ini h(Ti.^[T,x,Dx(p{T,x),u\ — u(T,x),e) > ( resp. <) 0. 



We have similar equivalent definitions of viscosity super and subsolution to (j4.ip in terms 
of test functions. 
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From the Markov property of the jump-diffusion process X, and uniqueness of a minimal 
solution Y to the BSDE (|230]) . we see that Yt = v{t,Xt), 0<t<T, where 



v{t,x) := , {t,x) £[0,T]xW, (4.4) 

is a deterministic function of (t, x), {Xl'^ , t < s < T} is the solution to (|2.1|) starting from 
X at time t, and {Ys*'^,i < s < T} is the minimal solution to ()2.10p - ()2.1ip with Xs = Xt"^ , 
t < s <T. Similarly, we define the function 

Vnit,x) := y^"'*'", (t,rE) G [0,T] xM^ (4.5) 

where {(y/*'*'^, Z"'*'^, C/"'*'^(.)), t < s < T} is the unique solution to (|3.ip with = Xt'^ , 
t < s < T. The relation between the penalized BSDE (j3.ip and the penalized IPDE (j4.2p 
is well-known from the results of 0. Although our framework does not fit exactly into the 
one of [2] , by mimicking closely the arguments in this paper and using comparison theorem 
in |20j . we obtain the following result. 



Proposition 4.1 The function Vn in (|4.5p is a continuous viscosity solution to (|3.ip . 

We slightly strengthen Assumption (HI) or (H2) by 

(HI') There exists a quadruple {Y, Z, K) e S"^ x L'^{W) x satisfying (IXT^ . with Yt = 
v{t,Xt), < t < T, for some function deterministic v satisfying a linear growth condition 

\v(t, x)\ 

sup ' ^ ^' < +CX3 

(t,x)G[0,T]xR<* 1 + FI 

(H2') There exists a quadruple {Y,Z,K,U) e x L^iW) x L'^ifl) x satisfying 
(|2.10p - ()2.1ip . with Yt = v{t,Xt), < t < T, for some function deterministic v satisfying a 
linear growth condition 

\v(t,x)\ 

sup V"^rT < 

Remark 4.2 Assumption (H2') (or (HI') which is weaker than (H2') in the case h{u,e) 
= —u) ensures that the function v in (|4.4p satisfies a linear growth condition, and is in 
particular locally bounded. Indeed, from (|3.1Up and by passing to the limit by Fatou's 
lemma for v(t,x) = Y^'^ = limY^"'*'^, we have 



t 

T 



sup |7;(t,x)p < C7(i + E|5(X^'^)|2 + e[ /" + E[ sup |i}(s,X*'^i2 



te[o,T] ^ ^Jt 

The result follows from the standard estimate 



E[ sup |A:*'^P] < C(l + |x|2), 

t<s<T 

and the linear growth conditions on g and v. 

By adapting stability arguments for viscosity solutions to our context, we now prove 
the viscosity property of the function v to (j4.ip . 
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Theorem 4.1 Under (H2') (or (HI') in the case : h{u,e) = —u), the function v in ()4.4p 
is a viscosity solution to (j4.ip . 

Proof. From the results of the previous section, we know that v is the pointwise hmit of 
the nondecreasing sequence of functions {vn)- By continuity of Vn-, we then have (see e.g. 

m p- 91) ■■ 

V = = hm mi^Vm where hm inf^,7;„(t, x) := hminf Vn{t',x'), (4.6) 

n-^oo n.^oo n ^ oo 

t — > t, X — > X 

V* = lim snp^Vn: where lim snp^Vnit^x) := limsup Vn{t' ^x'). (4.7) 

71 — * OO 
t — ' t, X — > X 

(i) We first show the viscosity supersolution property for v = v^. Let {t,x) a point in 
[0, T) X W^, and (p, q, M) G J~v{t, x). By ()4.6p and Lemma 6.1 in [7], there exists sequences 

nj oo, {pj,qj,Mj) G J~Vn,{tj,Xj), 

such that 

{tj,Xj,Vnj{tj,Xj),pj,qj,Mj) {t,x,v{t,x),p,q,M). (4.8) 

We also have by definition oi v = v^: and continuity of 7 : 

f (t, a; + 7(x, e)) < limmfvnAtj,Xj + j{xj,e)), \/eGE. (4.9) 

j-+oo 

Moreover, from the viscosity supersolution property for Vnj , we have for all j 

-Pj - - ^tr(cro-^(xj)Mj) - f{xj,Vnj{tj,Xj),a''{xj)qj) 

-nj h~{n''[tj,Xj,qj,Vn^] - Vn^{tj,Xj),e)X{de) > 0. (4.10) 
Je 

Let us check that the following inequality holds : 

mi h{n''[t,x,q,v]-v{t,x),e) > 0. (4.11) 
We argue by contradiction, and assume there exists some cq £ E s.t. 

h{v{t,x + 'y{x,eo)) + c{x,v{t,x),a'' {x)q,eo) - v{t,x),eo) < 0. 



Then, by continuity of o", /i, 7, c in all their variables, (j4.8p . (j4.9p . and the nonincreasing 
property of h, one may find some e > and some open neighborhood Oq of cq such that 
for all j large enough : 

h{vnj{tj,Xj +-f{xj,e)) + c{xj,Vn,i{tj,Xj),a''{xj)qj,e) - Vnj{tj,Xj),e) < -e, Ve G Oq. 
Since the support of A is E, this implies 

h~{n''{tj,Xj,qj,Vnj) - Vnj{tj,Xj),e)X{de) > eA(C'o) > 0. 
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By sending j to infinity into (|4.10p . we get the required contradiction. On the other hand, 
by (|4.10p . we have 

-Pj - {K^j),Qj) - ^tr(cro-^(xj)M,) - f{xj,Vnjitj,Xj),a''{xj)qj) > 0, 

so that by sending j to infinity: 

-p- {b{x),q) -^tr{aa^x)M)- f{x,v{t,x),a'{x)q) > 0, 

which proves, together with (|4.1ip . that f is a viscosity supersolution to (j4.ip . 

(ii) We conclude by showing the viscosity subsolution property for v*. Let (t, x) a point in 

[0,r) X R'^, and {p,q,M) e J+v*{t,x) such that 

M h(n^[t,x,q,v*]-v*(t,x),e) > 0. (4.12) 

eG-B 

From (j4.7p and Lemma 6.1 in 0, there exists sequences 

rij oo, {pj,qj,Mj) G J+Vnj{tj,Xj), 

such that 

{tj,Xj,Vnj{tj,Xj),pj,qj,Mj) {t,x,v*{t,x),p,q,M). (4.13) 

By continuity of the functions c, 7, and definition of v* , we also have 

limsu];)n''[tj,Xj,qj,Vnj] < Ti'^it, x, q,v*], Me£E. (4.14) 

j-+oo 

Now, from the viscosity subsolution property for v^^ , we have for all j 

-Pj - {H^j),Qj) - ^tr(crc7^(xj)M,) - f{xj,Vnj{tj,Xj),a''{xj)qj) 



-Uj h {n''[tj,Xj,qj,VnA - VnAtj,Xj),e)X{de) < 0. (4.15) 
Je 

Prom (j4.12p - (j4.13p - (j4.14p . continuity assumptions on h,c, and the nonincreasing property 
oi h, we have for j large enough 

h{7i''[tj,Xj,qj,Vnj] - Vn^ itj,Xj), e) > 0, Ve G E, 

and so 

/ h~{n''[tj,Xj,qj,VnA - Vnj{tj,Xj),e)X{de) = 0. 
Je 

Hence, by taking the limit as j goes to infinity, into (|4.15p . we conclude that 
-p- {b{x),q) - ^tv{aa'{x)M) - f{x,v*it,x),a'{x)q) < 0, 

which shows the viscosity subsolution property for v* to ()4.ip . □ 
We next turn to the boundary condition. 
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Theorem 4.2 Under (H2') (or (HI') in the case : h{u,e) = —u), the function v in ()4.4p 
is a viscosity solution to (j4.3p . 



In order to deal with the possible jump at the terminal condition, we need the following 
dynamic programming caracterization of the minimal solution. 

Lemma 4.1 Let {t,x) G [0, T) x R'^, and (Y^<^ ^ z^'^ ,U^'^ , K^'^) be a minimal solution to 
(j2.10p - (j2.1ip on [t,T] with Xg = Xl'^ . Then for any stopping time 6 valued in [t,T], 

is a minimal solution to : 

Ys = v{9,Xl'^)+ f f{X'^^-,Yr,Zr)dr + Kl'''-Kl^- (4.16) 

J s 

- j\zr,dWr) - f {Ur{e) - c{xl'^ ,Yr- , Zr, e)) ^l{dr, de) 



with 



h{Usie),e)>0 dF (g) dt X{de) a.e. on n x [t,e] x E. (4.17) 



Proof. Let Y^ be the minimal solution on [t, 9] of (I4.16p - (l4.17p (the existence of a minimal 
solution in the case of a random terminal time is obtained by similar arguments to those 
used in the case of a deterministic terminal time). For each u £ il., there exists a minimal 
solution y2,<^ on [9{u}),T] to ([2J0]) - (pnTT) . We then have from the Markov property of X 
that Yg^^-^ = V (^9(uj), Xg^^^{uj)^ for all uj £ Q. By a measurable selection theorem, there 

exists G such that for all lj G we have Ygj^^-^{uj) = Yg(^-^ = v (^9{uj), Xg^^-^{uj)^ 

and Yg{Lu) = Ys''^{lo) for s G [9{uj),T]. We then define the process Y by y|[t,0] = Y^ and 
Y\^e,T] = Y^. Hence, F is a solution on [t,T] to (l2J0]) - (l2lT] ). which implies Y > y*'^. 
Moreover, since y^'^ = v{9,Xg'^), it follows that (Yg'^ , Zl'^ ,Us'^ , Ki'^)g^[i^gj is a solution 
on [t,9] to KT6h - iKTl} i. Hence < y*'^' on [t,9], and therefore y^ = y*'^' on [t,9]. □ 

Proof of Theorem 14.21 (i) We first prove the supersolution property of u^, to (j4.3p . Let 

X G W^, and {p,q,M) G J~v^{T,x). By same arguments as in (|4.1ip . we have 

ini h{n^[T,x,q,v^]-v^{T,x),e) > 0. (4.18) 

Moreover, since the sequence of continuous functions {vn)n in nondecreasing and Vn{T, .) = 
(7, we deduce that w*(T, ■) > g, which combined with (|4.18p . proves the viscosity supersolu- 
tion property for to (j4.3p . 

(ii) We next prove the subsolution property of v* to (|4.3p . We argue by contradiction and 
asume that there exist xq gR"^, ip £ C^'2([0, T] x M") such that 

= (v* -ip){T,xo) = max {v* - if) (4.19) 

[0,T]xR'' 

and 

=: 2e > 0. 



mm 



ip(T,xo) - g{xo) , inf /i(W^[r,xo,-Di(/:'(T,xo),w*] - </?(r,xo),e) 
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By the upper semicontinuity of v* , the continuity of (p and its derivative, and the nonin- 
creasing property of h, there exists an open neighbohood O of (T,xq) in [0, T] x W^, and 
A,r > such that for all {t, x, a, /3) G O x {—A, A) x B{0, r), we have 



e < min 



f{t,x) - a- g{x) , 



(4.20) 



inf h{v*{t, X + ^{x, e)) + c(x, (/^(i, x) — a, a'' {x)[Dxip{t, x) + f3]) — [^{t, x) — a], e) 



eS-E 



Let {tk,Xk)k be a sequence in [0, T) x M'^ such that 

{tk,Xk) ^ {T,xq) and ■u(tfc, x^) ^ 7;*(r, xq). 



(4.21) 



Fix then 5 > such that for k large enough: [tfe,T] x B{xk,S) C O, and let us define the 
functions ipk by 

|2 



ipk[t,x) = ip{t,x) + C p KL-fc' 



+ VT-t, 



where < C < ^ A (5r, e 
oo, and Ck > is a constant to be chosen below. By (j4.19p . we notice that 



satisfies 4'\b{o,i) = 0, <AIb(o,i)= > ^^"^ 1™ 



x|-*oo i+ui 



(■;;* -(/jfc)(t,x) < -C ioT {t,x) £[tk,T]xdB{xk,6), 
and from the conditions on (p, we can choose Ck (large enough) so that 

c 



{v* - (fk){t,x) < 



for (t,x) G S(xfc,<5)^ X [t;^,,r]. 



(4.22) 



Since §-Ay/T^t) — > — c« as t X r, we have for /c large enough : 



d^k 
' dt 



Cipk{t,x) - f{x,(pk{t,x) - a,a''{x)D:c(Pk{t,x)) 



> for (t,x,a) G [tfc,r) X S(xfc,5) X (-A + C,A). 



(4.23) 



Fix now a* G (0, ^ A | A e), and let us denote = inf {s > ; / ^s-}^ = 
inf {s > tfc ; Xg ^ B{xk,6)] ATkAT where X'^ = X^*"^'' . Let us then define the quadruples 
{Y\Z\U\K^) on [tk^Ok] by : 



U^e) = v*{s,X^^+j{X^^,e))+ciX^.,Ms,X^.)-a*,a^{X^.)D,Ms,X^-)) 
- [ipk{s,X^^)-a*], 



and 



|^(r,X,^) +/:(^fc(r,X,^) + fiX^,Mr,Xh - a* iX^)D,Mr,X^))^ dr 




{ipk - a* - v*){r, X^^ + jiX^^ , e))fiidr, de) 



+ {^k{ek,Xl) - a* - v{dk,Xl)) 



31 



By construction and from Ito's formula on (pk{s, X^), we see that {Y^, Z^, U^, K^) satisfies 
(|4.16p on \tk,Of^]. From (j4.20|) . it is clear that the process satisfies the constraint : 

h{U^{e),e) > 0, dP ^ dt ^ X{de) a.e. on n x [tk,0k] x E. 

Observe also that 

M6k,X^J-a* > viOk^Xl) (4.24) 

Indeed, we have two cases: 

• Ok <T: in this case {Ok-, X^J ^ O, and since a* < |, we have by (|4.22p . 

iPk{ek,x^j-a* > v*{ek,x^j > viOk^xl). 

• ek = T: in this case {Ok,X^J = {T,X!^) G O. Since a* < e, we have by (|i:20]) 

M0k,X^J-a* > ^{9k,X^J-e>g{X^) = v{ek,X^J. 

Let us then check that K'' is nondecreasing on [t^, 9^]- First, on [t^, 0^), we notice that K'^ 
consists only in the Lebesgue term dr, and so is nondecreasing by (j4.23p . Moreover, we see 
that Kn > . Indeed, there are two possible cases: 

• ek<Tk: then K^^ = + M0k,X^J -a*- viOk^X^^), and by (gSll), we have 

• ek = Tk: then = - {^k{Ok,X^^) -a*- v*{ek,X^J) + {MOk,X^J - a* - 
v{ek,Xl)), and soKl > 

Therefore, the quadruple Z'"' , K^) is a solution on [tk.Ok] to (!4.16p - (|4.17p . and by 
Lemma |4.H we deduce that for all /c, 

ipk{tk,Xk) - a* = ip{tk,Xk) + \/T - tk - a* > v{tk,Xk). 
We finally obtain a contradiction by sending to cxo. □ 

4.2 Uniqueness result 

This paragraph is devoted to a uniqueness result for the QVI (j4.ip - (j4.3p . We need to impose 
some additional assumptions. 

(H3) There exists a nonnegative function A G ^^(M'^) and a positive constant p satisfying 

(i) CK + f{.,k,a^Dk) < pA, 

(ii) MeeE HT-CAix) - A(x), e) > for all x G R'^, 
(in) A(x) > g{x) for all x G M"', 

(iv) lim|^.|_^ TTR = oo- 
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Assumption (H3) essentially ensures the existence of a suitable strict supersolution 
to (j4.ip . We shall give in paragraph [6] some sufficient conditions for (H3). This strict 
supersolution allows to control the nonlocal term in QVI (j4.ip - (j4.3p via some convex small 
perturbation. Thus, to deal with the dependence of /, c on y, z, we also require some 
convexity conditions. 

(H4) 

(i) The function /(x, ., .) is convex in {y, z) £RxW^ for all x G M'^. 

(ii) The function h(., e) is concave in u G M a for all e (z E. 

(iii) The function c(x, ., ., e) is convex in {y, z) G M X for all (x, e) G M'^ x E. 

(iv) The function c{x, ., -z, e) is decreasing in y G M for all {x, z, e) G M*^ x M'^ x E. 

Theorem 4.3 Assume that (H3) and (H4) hold, and let U (resp. V ) he a Isc (resp. 
use) viscosity supersolution (resp. subsolution) to M.ip -f jlTSj) satisfying a linear growth 
condition : 

\U{t,x)\+\V{t,x)\ vy.^rnTl 
sup i — i < oo, Vi G [U, J J. 

Then, U >V on [0,T] x R'^. Consequently, under (H2') (or (HI') in the case : h{u,e) = 
—u), (H3) and (H4), the function v in (|4.4p is the unique viscosity solution to (j4.ip - ()4.3p 
satisfying a linear growth condition, and v is continuous on [0,T) x W^. 

Proof. • Comparison principle. As usual, we shall argue by contradiction by assuming 
that 

sup {V -U)> 0. (4.25) 

[0,T] xM'' 

1. For some A > to be chosen below, let 

U{t,x) = e^P^^'^*U{t,x) , V{t,x) = e^P+^'^'Vit^x) and ~k{t,x) = e^P+^'^' A{x) . 

A straightforward derivation shows that U (resp. V) is a viscosity supersolution (resp. 
subsolution) to 



mm 



dw 

pw - — - Cw - f {■,w,a'' D^w) , (4.26) 



inf h ■,7{'^w — w,e 

eG-B 



0, on [0, T) X 



t>d 



mm 



w{T^,-)-g , inf h{T,n^uj{T-,-)-w{T-,-),e) = on M'^ (4.27) 

eG-E 



where 

f{t,x,r,q) = e(^+^)*/(a^,?^e-(^+^)*,ge-(''+^)*) - Ar 
h{t,r,e) = e(''+^)*/i(e-(^+^)V,e), g{x) = e''P+^'>^ g{x) 
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and 

Hw{t,x) = w{t,x + ^{x,e)) + c{x,'w{t,x),a'' {x)Dxw{t,x),e) 

with 

c{t,x,r,q,e) = e^^+^^^^x, e'^/'+^^V, e"(''+^)*g, e) 

for all {t, X, r, q, e) G [0, T] x M'^ x M x R'^ x E. Since / is Lipschitz, we can choose A large 
enough so that / is nonincreasing in r. Denote W = (! — //)[/ + /iA with > 0. By (j4.25p 
and the growth condition (H3)(iv) of A, we have for /x small enough 

sup {V -W) = {V -W){tQ,XQ)>^. (4.28) 

[0,T]xR'* 

for some (^052:0) G [0,7"] x IR'^- Moreover from the viscosity supersolution property (|4.26|) - 
(|4.27|) of U , and the conditions (H3)(i), (ii), (H4)(i), (ii), (iii), we see that is a viscosity 
supersolution to 

pw - Cw - f {-.w^cr^D^w) > 0, on [0, T) X M"^, (4.29) 
ot 



inf h{-,n''w -w,e) > fiq, on [0, T] x M'^, (4.30) 
where q{t,x) = e^P+^'^HnleeE HWAix) - A{x),e) is positive on [0,r] x M"' by (H3)(ii). 

2. Denote for all [t, x, y) G [0, T] x R"' x and n > 1 

@n{t,x,y) = V{t,x) -W{t,y) - ipn{t,x,y), 

with 

ipn{t,x,y) = njx - + |x - xol"^ + |t - top- 
By the growth assumption on U and V and (H3)(iii), for all n, there exists {tn,Xn,yn) S 
[0, T] X M*^ X M*^ attaining the maximum of 6„ on [0, T] x M'^ x M'^. By standard arguments, 
we have : 

{tn,Xn,yn) ^ {to,Xo,Xo), (4.31) 

n\xn-yn\^ ^ 0, (4.32) 
)-W{tn-,yn) V{tQ,XQ) -W{tQ,XQ). (4.33) 

3. We now show that for n large enough 

inf /l(t„,?Y^[t„,X„,i:>a;V'n(tn,a;n,yn),t^] - t^(tn,a;n),e) > 0. (4.34) 

On the contrary, up to a subsequence, we would have for all n, 

\vdh{tn,H''[tn,Xn,D:^Lpn{tn:Xn,yn),V]-V{tn,Xn),e) < 0, 
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and so by uppersemicontinuity of V, compactness of E, there would exist a sequence (cn) 
in E such that 

h{tn,'H''"[tn,Xn,D^ipn{tn,Xn,yn),V]-V(tn,Xn),en) < 0. 

Moreover, by the viscosity supersolution property of W to ()4.30p . we have 

h{tn,Ti.''"[tn,yn,-Dyipn{tn,Xn,yn),W] -W{tn,yn),en) > fJ.q{tn,yn)- 

Prom the nonincreasing and the Lipschitz property of h{.,e), we deduce from the two 
previous inequahties that there exists a positive constant ij such that 

n J Xn , yn. ni y-n) 

<^ ?Y " \tfi , Xn , Dx (tn i Xn , yn ) ? ^] ^ (^n; ) ; 

which is rewritten as 

V{tn,Xn) - W{tn, yn) + n: yn) 

))-W{tn,yn + l{yn,en)) + ^Cn (4.35) 

where 

^Cn — C ^tn, XnjV (tnj Xn) T O'^ (^Xn)Dx'-Pn(tnT Xnj yn) T (^n^ 

- c(^tn,yn,W{tn,yn),-Cr''{yn)Dyipn(tn,Xn,yn)^ ■ 

Now, we write AC„ = AC^ + AC^, with 

C ^tn; ^nj ^^(tm yn)j {Xn) Dx^fnitrn ^n; 2/71)5 ) 
ACl = c(tn,Xn,Witn,yn),Cr''iXn)Dx(Pn{t n ) Xn 1 yn ) i I 

C ^tn; Xm {tni 2/n) j '^^ iyn)-Dy(pn (tm Xni yn)i 
AC^ = c{tn,Xn,W{tn,yn),-Cr'' {yn)DyLpn{tn,Xn,yn),ei?j 

C ^tnj J/nj U itn^ yn)) '^^ {yn)E)y^n (tn j j ^/n) ) • 

We have V{tn,Xn) - W{tn,yn) {V - W){to,xo) > by and (|05D . Hence, for 

n large enough, y(t„,x„) > VF(tn)yn)i and so from the nonincreasing condition (H4)(iv) 
of c, we have AC^ < 0. Since a'' {xn)Dx(pn{tn,Xn,yn) + a'' {yn)Dyipn(tn,Xn,yn) ^ by 
the Lipschitz condition on a and (j4.32p . we deduce with the Lipschitz condition on c that 
limsup„^oQ ACn < 0. By (j4.3ip and continuity of c, we have lim„^oo ^C'n = 0- Therefore, 
we obtain 

lim sup ACn < 0. 
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Up to a subsequence, we may assume that (e.„) converges to cq in E. Hence, by sending n 
to infinity into (j4.35p . it follows with (j4.33p and the upper (resp. lower)-semicontinmty of 
V (resp. W) that : 

{V -W)ito,xo + ^{xo,eo),xo + 'j{xo,eo)) > {V - W){to,xo) + m{to,xo) 

> {V-W){to,xo), 



a contradiction with (|4.28|) . 



4. Let us check that, up to a subsequence, tn < T for all n. On the contrary, tn = to = T 
for n large enough, and from (j4.34p . and the viscosity subsolution property of V to (|4.27p . 
we would get 

V{T,Xn) < g{Xn). 

On the other hand, by the viscosity supersolution property of U to (j4.27p and (H3)(iii), 
we have W{T,yn) > giUn), and so 

V{T,Xn)-W{T,yn) < g{xn)-g{yn). 

By sending n to infinity, and from continuity of g, this would imply {V — W){tQ,XQ) < 0, 
a contradiction with (|4.28p . 



5. We may then apply Ishii's lemma (see Theorem 8.3 in j7]) to (t.„, y„) G [0, T) 

that attains the maximum of 0„, for all n > 1 : there exist (p^, g^, M„) G J'^'^V{tn, x„) 
and {p'^,q'^,Nn) G P'-W{tn,yn) such that 

P^-P^ = dtiPn{tn,Xn,yn) = 2(t„-to), 

Qy — Dx'^ni.tm Xniyn) 1 Qy^ — Dyipnitn^ Xn, yn) ^ 

and 

where An = D'^^ y-^^nitn,Xn, yn)- From the viscosity supersolution property of W to (|4.29p . 
we have 

pW{tn,yn) -P"^- {b{yn), Dyip{tn,Xn,yn)) " ^tr (cr(y„)(T^ (y„,) A^n) 

-fitn,yn,W{tn,yn),-Cr''iyn)Dyip{tn,Xn,yn)) > 0. 

On the other hand, from (j4.34p and the viscosity subsolution property of V to ()4.26p . we 
have 

pV{tn,Xn) -P'^ + {b{Xn),Dr,ip{tn,Xn,yn)) " ^tr(cr(x„)cr^ (x„)M„) 

-f{tn,Xn,V{tn,Xn),(T''{Xn)Dxip{tn,Xn,yn)) < 0. 
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By subtracting the two previous inequalities, we obtain 

- {b{Xn),D^ipn{t )) + {b{yn),Dyipn{tn,Xn,yn)) 

+ ^tr {a{xn)a'{xn)Mn - a{yn)(T'{yn)Nn) , (4.37) 

where 

- f(tn, y-n, W{tn, y-n), -O"^ (y„)-D.y Xn,yn))- 

From (j4.3ip . we have ~ — > as n goes to infinity. From the Lipschitz property of b, 
and (|4.32p . we have 

lim({b{Xn),Dx(pn{tn,Xn,yn)) + {b{yn),DyVn{tn,Xn,yn)] = 0. 

As usual, from (|4.36p . (|4.3ip . (|4.32p . and the Lipschitz property of a, we have 
limsuptr ((t(x„)o-^(x„)M„ - (j(7/„)(T^(y„)iV„) < 0. 

Moreover, by the same arguments as for c, using the nonincreasing property of / in its 
third variable, and the Lipschitz property of /, we have 

limsupAi^„ < 0. 

n— >oo 

Therefre, by sending n — > cxo into ()4.37p . we conclude with ()4.33p that p{V — W){to,xo) < 
0, a contradiction with (j4.28p . 

• Uniqueness for v. The uniqueness result is then a direct consequence of the comparison 
principle, and the continuity of v on [0, T) x M.'^ follows from the fact that in this case = 

V*. □ 

Remark 4.3 As a byproduct of the comparison principle in Theorem 14.31 we get the 
continuity of the value function v on [0, T) x M.'^. Since the jump-diffusion process X is 
quasi-left continuous, then so is the minimal solution Yt = v{t,Xt) to the BSDE with 
constrained jumps, and the penalized approximation = Vn{t,Xt). This implies that the 
predictable projections ^Y and pY"^, respectively of Y and Y^, are equal to ^Yt = Yj- and 
pyn ^ yn _ Therefore, Y^ = hm 

n—^co • From the weak version of Dini's theorem, see [9] 
p. 202, this yields the uniform convergence of Y"' on [0, T], i.e. lim^^oo supjg^ \Y-P — Yt\ 
= 0, and so by the dominated convergence theorem, the convergence of Y^ to y in 5^ : 

lim ||y"-y|L2 = 0. (4.38) 

n— >oo 'S 

Then, by applying Ito's formula to — ^"P, and by using similar arguments as in 

Lemma 13.31 one can show that 



\^n+p _ ^nu2 _|_ II rr"+P _ rr" II 2 _(_ || _ ?^'"||2 < /^||y"+P _ y"||2 

I IIl2(w) " II II52 - II Ils2- 
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Together with (j4.38p . this imphes that (Z""),^ (U"')n i^"')n are Cauchy sequences 
respectively in the Banach spaces L^(W), L^(/2) and 5^. Therefore, under the additional 
conditions (H3) and (H4) with respect to Theorem 13. 11 we obtain the strong convergence 
of (Z", U^') in L^(W) X L^(/2). Notice also that in this case, the limiting process K of 
in 5^ is not only predictable but inherits the continuity property of K"^, see also Remark 

E31 



5 Numerical issues 

In this section, we formally discuss the numerical implications of our representation results 
and approximation by penalization for QVIs (|4.1|) . We first briefly recall the classical 
approach for numerically solving the QVI (II. ip arising from the impulse control problem 
(jl.4p . This is based on an approximation by iterated free boundary or optimal stopping 
problems : Starting from the function 

T 



uo{t,x) = Eb(X°'*'^') + ^ fix^/n], 



solution to the Cauchy problem : 

-^-Cuo-f = 0, on[0,T)xM'^, uo{T,.)=g on 

we construct the sequence of functions by induction as the solution to the optimal 
stopping problem : 

u„^,{t,x) = sup E[Wn„(T,xO'*'")], 



which satisfies the obstacle PDE 



mm 



dt 



^u„+i - f , Un+i - = 0, on [0,r) X M , n„+i(r, .) = Tiun on 



where Ti is the nonlocal operator defined in (jl.3p . We refer to the book [16] for a more 
detailed description of this approximation scheme. Such an numerical approach is com- 
putationally demanding, since it requires at each induction step n, the resolution of an 
optimal stopping problem. Moreover, at step n + 1, for determining the function u^^-^ at 
one point, one needs to compute the function u„ in the whole space due to the nonlocal 
term in the obstacle Tiun- 

We consider the general QVI (j4.ip and we propose here a numerical approach based on 
the probabilistic representation of the solution to this QVI by the constrained BSDE ()2.10p . 
We only describe the steps of the algorithm and postpone the analysis of the convergence 
to a future research. 

Step 1. Approximation by penalized BSDE. That is, we use (13. ip to approximate 
(iZTnll . The convergence of (Y", Z", /sT") to {Y,Z,U,K) is due to Lemma [M] and 
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TheoremO We note that, by denoting Vp{e) := U2{e)-c{Xs-,YJl, Z2_,e), (Y", Z", F") 
satisfies the following BSDE: 

rr = giXT) + r [ fniXs,Yp,Z'^,VPie),e)Xide)ds (5.1) 

Jt JE 



t JE 



y,"(e)A(de,ds), 



where 



fn{x, y, z, V, e) := J^fi^^ y,z) -v + nh {v + c{x, y, z, e), e). 



(5.2) 



Step 2. Discretization in E. For each m, let EY^,---,El^ be a partition of E the 
state space of the jump size. Denote ef := eX{de). Let (X"^, F'^''", Z'^'™, F"'"^) 

denote the solution to the following BSDE: 



Y 



nt ft ™ l-t 

x+ / h{X^)ds+ / a{X^)dWs + y] / 7(^r,ef)/x(ds,£^f ); 
JO Jo j-^^ ./o 

g{XT) + Y,KET) / /„(X„n"'-,Z,"'™,r^'-(j),e7)ds (5.3) 

.7 = 1 ''^ 



where 



t JE 



\{Ef) 



Vr'"'ie)fl{ds,de), 



Vr'"'ie)X{de). 



By the Lipschitz conditions, one can easily see that (j5.3p is well-posed. Moreover, one can 
show that once maxi<j<m diam(£'™') — > where diam(£'™) := sup{|ei — 62] : ei,e2 £ -^j"}) 
then 



lim 

m-->oo _ 



1-^ ~ Ils2 + I|j^ + l|Z -z IIl2(W) 



0. 



Step 3. Discretization in time. This is an extension of the work by Bouchard and Elie 
J5j, which studies the case m = 1. For any N, let h := jr and ti := ih,i = 0, ■ ■ ■ , N. First, 



define XV^'^ := x, and for z = 0, • • • , — 1, 



'-to 



Next, define Yj]]^'^ := g{X^^), and for i = — 1, • • • , 0, 



Z 



n.m,N 



rn,m,N / ■\ 



^i, =Etjyt^^^' ] + h2^x{Ej )fn{x^^' .y^/ ,r^; (j),ej. )(5-6) 

i=i 



/iA(£;: 



i-E,[y,:^'^/i((t„t,+i]xi?-) 



n.m,N-\ 



,J = 1 



! 5 ""5 



m; 



(5.4) 
(5.5) 
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Here Ef . denotes the conditional expectation under Tt^ ■ Notice that y^"'"^'^ ig defined via 
implicit scheme. One can also define it via explicit scheme by replacing the y^"'"^'^ inside 
fn with ^"l™' • Following the arguments in [5], one can prove the convergence of this 
time-discretization approximation : 



lim \ max E 

iV^oo I 0<i<A'-l 



+ 



{ 



Y'ft.lll, T 7-11 

ti <t<ti^l 



n,m, N i2 



N-1 



1=0 



ti + l 



n,m iyn,m,N \2 , 



Step 4. Approximation of the conditional expectations. The last step consists in 
the approximation of the conditional expectations arising in (j5.4|) - (|5.5p - (j5.6p . There are 
several approaches proposed in the literature. We adopt here the Longstaff-Schwarz [T3] 
method by least square projection and Monte Carlo simulations. By induction one can 
easily see that 



rn,m,N n,m,N / ^m,N \ yn,m,N in,m,N / ^m,N \ j^n,m,N , n,m,N / ^''TT-y^ 



n,m,N !n,m,N i n,m,N 



i/j-' ' . For any L, choose some basis 



for some deterministic functions v- 

functions {v^, 4>f, V'f ), I = - ■ ■ ,L. For any M, simulate M independent copies of ^ 
W,^: and fl'^iiU, ti+i] xEJ^),k = l,---,M. 

First, for each k, define X^^^'^^ := x, and for z = 0, • • • , A'" — 1, 



X 



m,NM 



m 



m,N,M \ 



k,t 



i+i 



Next, define Y^'^'^'^'^ := g{X^^^'^) for each k, and for i = N — 1, ■ ■ ■ ,0, we define 

/ r^n,m,N,L,M -r^n,m,N,L,M -t^n,m,N,L,M\ r ii 
y^k,U ' k,U ' ^k,U ) ^ tOilOWS. 



M L 

(«„... , a,) := arg^min ^ E | ^^^Xf' '"'[^^ " " E 
' ' k=l 1=1 



yn,m,N,L,M 



1=1 



(/3i(j), • • • , PlU)) ■= arg min — x 



X] \ h\(T^ 



Y 

m\ ti+i 



n,m,N,L,M ~k 



^n,m,N,L,M 
k.ti 



((t.,t.+i]xi?f)-EAV'f(xr:f'*-') 



1=1 



1=1 



(7ir--,7L) :=arg min 

71. ■■■.7l M 



M L 
t \-,^n,m,N,L,M \ ^ L/ vrn,N,M\ 



k=l 



1=1 
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h.ti 



k,U > 



1=1 



Finally we solve the following equation for ^ j. ' ' ' via Picard iteration: 

m 

6 Some sufficient conditions for (H2') and (H3) 

In this section, we provide various explicit conditions on the coefficients model, which ensure 
that the general assumptions (H2') and (H3) hold true. 



6.1 Existence of the solution to BSDE with jump constraint 

We first consider a case where we have upper bounds for the coefficients and h{u, e) = —u. 

Proposition 6.1 Suppose that h{u,e) = —u, and assume that there exist real constants 
Ci , C2 and 7] such that 

g{x) < Ci + {r],x) , c{x,y,z,e) + {r],j{x)) < and f{x,y,z) + {rj,b{x)) < C2, (6.7) 

for all {x, y,z,e) xRxW^ X E. Then (H2') holds true. 

Proof. Let us define a quadruple {Y, Z, K, U) by : Yt = Ci + C2{T - t) + {rj, Xt) for t < 
T, Yt = 9{Xt), Zt = a{Xt^).v, Ut{e) = and 

Kt = {c2-vb{Xs)- f{X,,Ys,Zs)}ds 

- ^ J^{c{Xs-,Y-,Zs,e) + {v,7{Xs-))} fi{ds,de), t<T, 

Kt = Kt- + Ci + {7],Xt) -giXr). 

From (|6.7p . the process K is clearly nondecreasing. Moreover, from the dynamics of X, 
and by construction, we see that the quadruple {Y,Z,K,U) satisfies ()2.10p - ()2.13p and the 
function v{t, x) = Ci + C2{T — t) + rj.x clearly satisfies a linear growth condition. □ 

We next give an example inspired by [1] where the jumps of X vanish as X goes out of 
a ball centered in zero in the case of impulse control. 



Proposition 6.2 Suppose that h{u, e) = —u, f, c does not depend on y, z, and assume that 
c < 0, 7 = on {x G M"' : \x\>Ci}xE for some Ci > 0. Then, (H2') holds true. 



Proof. We consider the function v : 

vit,x) = supE'^Ux^")+ / + 



c{X^'^ , e)ix{ds, de) 
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Since c < 0, and the choice of u = 1 corresponds to the probabihty measure = P, we 
see that v < v < v where 



v{t, x) 
v{t, x) 



E 



c{X 'J^, e)fi{ds, de) 



sup 



g{X'/) + I f{Xr)ds 



t,x\ 



The function v clearly satisfies a linear growth condition by the linear growth conditions 
on f, c and the standard estimate for X. Moreover, under the assumptions on the jump 
coefficient 7, it is shown in [3] that v satisfies a linear growth condition. Therefore, v also 
satisfies a linear growth condition. 

Let us now define the process Yt = v{t,Xt), which is then equal to 



Yt 



ess sup E'' 



9{Xt) + / f[Xs)ds + 



T 



c{Xg- , e)ii[ds, de) 



Tt 



and lies in 5^ from the linear growth condition, ad the estimate (j2.2p for X. From Theorem 
we then know that there exists (Z, [/, K) £ L^{W) x L2(/1) x such that (Y, Z, U, K) 
is the minimal solution to (j2.10p - (j2.13|) . and so (H2') is satisfied. □ 

We finally consider a case for general constraint function h. 

Proposition 6.3 Assume that there exists a Lipschitz function w G C^(M'^) satisfying a 
linear growth condition, supersolution to (j4.3p . and such that 



{b,Dw) + ^tr{aa''D'^w) + f{-,w,a^Dw) < C, 



on 



for some constant C. Then (H2') holds true. 
Proof. Let us define a quadruple {Y, Z, U, K) by 

Yt = w{Xt) + C{T -t), t<T, YT = g{XT), 
Zt = a^{Xt-)Dw{Xt-), Ut{e)=w{Xt-+j{Xt-,e))+c{Xt-,Yt-,Zt,e)-w{Xt-), and 

Kt = / [C - {b{Xs),Dw{Xs)) - -tv{a{Xs)a^Xs)DMXs)} - f(.Xs,Ys,Zs)]ds, t < T, 
Jo ^ 

Kt = Kt~ + w{Xt) - giXr). 

From the conditions on w, we see that (Y , Z , K ,U) lies in x L^(W) x L^(/i) x A^. 
Moreover, by Ito's formula to w{Xt) and the supersolution property of w to (j4.3p . we 
conclude that {Y,Z,K,U) is solution to (l230D - (l2lT] ). and v{t,x) = w{t,x) + C(T - t) 
satisfies a linear growth condition. □ 
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6.2 The strict supersolution condition (H3) 

We give a sufficient condition for (H3) in the usual case where / and c do not depend 
neither on y nor on z. 

Proposition 6.4 Consider the case where h is given by 

h{u, e) = —u. 

Assume that there exists a constant a > such that 

-a < |a; + 7(j;,e)p - |xp V(x,e)eR'^xE 
P := inf , |o >0 

Then assumption (H3) holds true. 

Proof. We set A{x) := + C with C large enough so that A > g, i.e. (H3)(iii) is 

satisfied. A straightforward computation shows that 

inf hiWKix) - Mx), e) > a/3 > 

and hence (H3) (ii) is satisfied. Clearly, (H3) (iv) holds as well. Finally, it follows from the 
linear growth assumption on h and cr that (H3) (i) holds for a sufficiently large parameter 
p. □ 
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